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ABSTRACT 


The need for a proper theoretical explanation of the residual 
film, reported Ly Tarious workers, has led to the development of 
microcontinuum theory of lubrication* Even after a decade this 
theory seems to be in a preliminary state. Scores of lubrication 
problems are yet to be analysed from this view point. For example, 
little work has appeared on dynamically loaded bearings and porous 
bearings from the microoontinuum view point. Moreover, no attempt 
has so far been made to study the impact of the micropolarity 
on finite dimensional bearings and bearings with rough surfaces. 

Thus, in this thesis, an endeavour is made to study these 
lubrication problems using the micropolar fluid theory. The thesis 
is divided into IX Chapters, 

Chapter I is of introductory type, wherein are discussed 
important experimental evidences, which establish the existence of 
the residual film (rheological abnormalities) . It also emphasizes 
the fact that the origin of the so called rheological abnormalities 
may lie in the micro-mechanical motion of the fluid element. Thus 
this Chapter lays the foundation for the present worjt. 

In Chapter II, the basic eq.uations for the micropolar fluid 
are given. Attention is restricted to a three dimensional laminar 
flow of an incompressible micropolar fluid in a finite slider type 
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bearing, Application of usual lubrication assumptions and a subsequent 
detailed order of magnitude analysis reduces the governing equations 
to two systems of coupled ordinary differential equations to be solved 
for velocity distributions and microrotation velocities (no pa.per, 
published so faf, gives such a detailed order of magnitude analysis 
for the three dimensional problan) . Generalized Reynolds equation 
is then obtained in a conventional manner by integrating the continuity 
equation. Also given in this Chapter are expressions for flow fluxes 
and shear stresses , 

In Chapter III, the generalised Reynolds equation, obtained 
in Chapter II, is applied to analyse the problem of two dimensional 
squeeze film of a ball in a spherical seat and to some three dimensional 
r.on-cyclic squeeze films, in an effort to study the effects of rigid 
particle additives for the three dimensional micropolarity model. 

In Chapter IV, dynamically loaded short journal bearings are 
studied* Reynolds equation for the general cas», xvhoro both wedge 
and squeeze films are active, is derived. Detailed consideration is 
given to the dynamic behaviour of squeeze film in a short journal 
bearing subject to a sinusoidal loading with no journal rotation. 

Chapter V, also gives an analysis of dynamically loaded hearings* 
The problem considered is that of the roller bearing in combined rolling, 
sliding and normal motion, under cavitation boundary conditions. 



In Chapter TL, the micropolar theory is applied to the analysis 
of hearings v/ith rough surfaces. Generalized forms of Beynolds 
equation are derived for three cases of roughness, namely, longitudinal, 
transverse and isotropic, using the stochastic approach. These 
equations are subsequently applied to the problem of infinitely 
long journal bearing whose stationary surface (bearing) is assumed 
to be rough, in an order to study the interaction of micropolarity 
with surface asperities, using the half Sommerfeld boundary conditions. 

In Chapter YII, the generalized Reynolds equations, obtained 
in Chapter YI, are applied to study the squeezing between parallel 
plates. Three geometries, namely, infinitely long parallel plates, 
circular plates and finite rectangular plates are considered. 

In Chapter YIII, the problem of synovial joints lubrication is 
studied, by approximating the joints by a spherical bearing, the 
cartilage surface as porous (but rigid) and the lubricating synovial 
fluid as micropolar, A modified form of Reynolds equation is derived 
using the slip boundary conditions. 

Chapter IX is of conclusive nature wherein are presented the 
main conclusions of the work done in the thesis and the conclusions 
are oompareo with the various experimental evidences. The one overall 

conclusion that emerges is the existence of increased effective 

\ 

viscosity in the proximity of a solid surface. This increase is further 
enhanced when the surfaces are rou^. 
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CHiPTER I 


GEMERAL INTRODUCTION 


1.1 WHY LUBRICATION 

Lubrication plays a most vital role in our great and complex 
civilization. To estimate the importance of this role one need only 
consider that evexy moving part of every machine is subjected to 
friction and wear. Tliroughout the centuries from ox-ca3rfc axles to 
grinding wheel spindles, one of man’s most persistent problems has 
centered around reduction and control of friction and wear. Friction 
consumes and wastes energy. Wear causes changes in dimensions and 
eventual breakdown of the machine element and the entire machine and 
all that depends upon it. 

Today the constituent parts of modern high speed engines are 
frequently subjected to a very high rate of 7/ear, which can lead to 
serious operational disturbances, if not to the total deterioration 
and destruction of engines. The need for continuous engine operation 
today necessitates the timely removal and replacenent of worn engine 
parts, which can be a costly operation, A nation’s economy can be 
seriously impaired through excessive losses of engine and machine 
parts due to wear, 'so that an intensive campaign against friction and 
wear is in the best of interest of increased production and national 
efficiency. 

So whether one is concerned with reducing the waste of world’s 
production of energy, conserving a nation’s critical nattiral resoxirces, 



maintaining production schedules in a plant, creating a larger 
margin of profit in plant operation, or even mdcing the fa:nily 
automobile run longer and better, one must be interested in the 
study of lubrication. 

1.2 EARLf DAYS Ih LUBRICATIOl 

Lubrication, in one sense, is as old in civilization as the 
wheel and axles. When a tomb vras opened in Egypt, some years ago, 
one of the chariots still had soma of the original lubricant on the 
axle. This was analyzed and found to be sticky and slightly 
greasy [ 1 ] . It contained road dirt such as quartz sand, compounds 
of aluminum, iron, and lime. It had a melting point of 120'^F', which 
suggests that it might have been mutton or beef tallow, either of vj'hioh 
would have proved suitable for axle lubrication in that warm country. 

Until comparatively recent years all liibricants were largely 
of animal, vegetable, or marine origin : mutton tallow, lard, gooze 
grease, fish oils, castor oil, oottonceed oil and other vegetable 

oils, etc. The use of mineral oils dates back principally a little < 

1 

over a hundred years . i 

( 

Pioneering work in the field of lubrication dates back to the I 

year 1847» when -^dams [ 2 ] developed and patented several rather good | 
designs for railway axle bearings, liany of the plans shov/ a remarkable | 
compr^ension of the fundamentals of proper bearing design. 

The new science of lute icat ion, however, based on hydrodynamics 1 
and oil flow, was not recognized as such until 1883» when Petrov [ 5 J 
made the first significant attempt to analyze theoretically the friction 



effect of film lubrication. It nas then recognized that the film rather 
than the bearing material could be a prime consideration. 

In the same year that Petrov published his theory, unexpected 
experimental results were reported by Tower [ 4 ] . He demonstrated 
that a loaded, oil-lubricated journal bearing is subject to local 
pressures substantially hi^er than its mean pressure, Tliree years 
later Reynolds [5] » unaware of Petrov’s theory^ was able to explain 
Tower’s results. Prom the properties of a thin laminar fluid film 
he derived and published not only the descriptive differential eq.uation, 
which even today bears his name, but also certain solutions to this 
equation that agree well with the experimental measurements of 
Tower [4 ] • Several yeetrs later in 1904> Sommerfeld [6] , provided 
elegant theoretical extensions to the journal bearing lubrication 
problem • 

Of the thousands of papers on film lubrication which have 
appeared since 1886, in the wake of Reynolds olassioal work on lubrication 
theory, few have offered modification of the theory. In general, papers 
on film lubrication have extended the application of the fundamentals, 
as set forth by Reynolds, to a variety of bearing shapes with steady 
and transient films of variable viscosity, density and temperature* 

In 1928 , Boswall [ 7 ] summarized the theory of film lubrication to 
that date. By now, however, much additional information about lubricating 
films has been developed. A proper historical review of '^^’6 same 
would be a Herculean task and beyond the scope of the present thesis. 
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Inspite of the rapid growth of the lubrication theory in the 
"fcii 

20 century, the picture is far from complete. A perfect correlation 
between theory and practice has not yet been obtained. This is fundamental! 
because the basic assumptions frequently circumscribe a rather tight ciixsle 
within which the applicability of the theory is Talid, 

1*5 CLASSIFICATIOir OF LUBRICATIOIT 

In surveying the entire field of lubrication with its many 
vital physical and chemical manifestations, the one phenomenon that 
occurs most frequently and persistently and which probably has by far the 
most value is the creation of a load-carrying fluid film. For the builders 
and users of machinery and mechanism of various kinds, development of 
this fluid film is quite a natural action and serves them obediently 
and tirelessly if given even the slightest encouragement , It is obvious 
that much of our industrial life would grind to a screechisog halt without 
this lubrication phenomenon. 

The thickness of the fluid film not only determines between 
the various kinds of lubrication but also the regime of applicability 
of the theory. At one extreme is the thick film lubrication, where 
the bearing surfaces are completely separated by a continuous fluid 
film of the order of 10 cm. and' the bearing characteristics mainly 
depend on the viscosity of the lubricant in bulk. At the other extreme 
is the boundary lubrication, where a fluid film only a few molecular 
layer thick is left between the bearing surfaces. Between these two 
extremes, lies, the region of thin film or mixed lubrication. The fluid 
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film is very thin, of the order of 10“^ cm, so that the load is partly 
supported hy the roughness asperities and partly hy the fluid film. The 
region of thin film lubrication is an indeterminate state in which the 
strict laws of hydrodjmamic lubrication are no longer valid. Some of the 
most exciting research is being done in this field. It is actually 
a meeting place for a chemist and an engineer as v/ell as a metallurgist 
and a mathematician. It is in this region that the rheology plays an 
important role. It is in this region that the proximity of a solid 
surface drastically modifies the properties of the lubricant . It is 
also an area of great endeavour, where the experimental facts are to be 
correlated with theory. 

1.4. LUBRICATIOIir AMD RHEOLOGY 

Interaction between the fields of lubrication and rheology has 
be.cotne increasingly important within the past few years with the growing 
awareness of the complexity of the lubricant reaction to severe and ; 

transient stress conditions, the continued demand for well lubricated, 
high speed, highly loaded machine elements, extensive use of synthetic | 
lubricants and use of additives. | 

I 

The general definition of rheology is the science that deals ! 

t 

with the deformation and flow of materials. Similarly, one may state, ; 

that the primary objective of lubrication is to separate rubbing surfaces ' 

1 

! 

by a layer of lubricant to prevent excessive friction and wear. It is cleg 
that the lubricant during the lubrication process will be subjected to a 
quite high shearing stresses and velocity gradients. The rheological | 
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properties of the Ixihricants are, therefore, lihely to have an important 
effect on its action. 

In early lubrication uork, only the ITewtonie.n viscosity of the 
lubricant was important. However the applied stresses have increased 
and th-e range of environmentation conditions, under which the lubricant 
should work, have become more extreme • The fact that many mineral 
base oils also manifest non-lewtonian behaviotir under the severs shear 
conditions they often encotmter in practice also support this point 
of view , 

Non-Newtonian behaviour is also almost invariably observed in 
various lubrication systems, when additives are used. These additives 
are put into the lubricant for a variety of purposes and do a great ^ 

deal to improve the lubricating oils which nature and the refiner produce* 
Modern lubricating-oil additives, based on years of scientific research, 
designed to meet the extreme demands of modern machines and for high 
performance ratings, under actual working conditions, have become : 

indispensable in many applications. ! 

The ever growing trend to transmit greater power through ^ 

mechanisms smaller in size and weight imposes an increasing burden | 

on the lubricant. New combustion problems, greater surface loadings, ’ 

1 

wider range of operating temperatures and often greater stalling or sliding! 

i 

! 

speeds in bearings, all subject the lubricants to higher operating j 

I 

temperat\ires, loads and overall abuses, ' , 

i 

These examples show just a few areas where a knowledge of the 
rheological properties of the lubricant is lilcely to be important, Nor 
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a. proper comprehension, one must look for the experimaatal evidences 

been 

where such proiDerties have been found, or have/soufrht for and found 
lacking. Hence in the next section a "brief stmmary of such experimental 
studies is given. 

1*5 EKPEdTMEimL EyiPEIlCES 

The manner in which a liquid in contact vath a solid surface 
is influenced "by the proximity of that surface has "been the subject 
of intensive study and a matter of controversy. It appears obvious 
that the molecules of the liquid in the immediate vicinity or those 
adhering to the surface, or of the adsorbed layers on the siirface 
must necessarily behave different from the adjacent molecules of 
the liquid in "bulk. However the nature of this surface influence 
and the depth to which it penetrates is not yet properly understood 
and the experimental results have given rise to widely divergent 
opinions. It is generetlly accepted that the surface zone of a 
liquid is not merely a monomolecular layer but a mltiraolecular ; 

layer [8,9 ] » i.e. the molecular orientation extends effectively 

I 

to many molecular Imgths, I 

[ 

Tb:om a practical point of view, however, the important question i 

i 

is not whether deep molecular orientation exists, but whether it gives | 

I 

rise to rheological abnormalities in the neighbourhood of a solid surfacl 
i.e, - whether a fluid in the immediate vicinity of a solid surface, devej 
properties significantly different from the adjacent molecules in bulk. 
The experimental evidences are viewed from this angle, I 
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KingsTaurj [ 10 j seems to be a,mong the earliest investigators 
of the question of existence of r’.ieologioo,! abnormalities in proximity 
of a solid surface. He observed an enhanced viscosity in that part of the 
fluid which is in the region of attraction of the surface molecules of the 
metals. This enhancement may range from a small increase in viscosity 

to rigidity and the depth of the affected zone ranges from 10 mm to 

“5 

10 ram. 

Hardy and Hottage [ 11-15 ] offer proofs that the liquid film 
must at least be of 0,007 'to 0,010 mm in thickness, if any of the liquid 
is to be beyond the range of surface influence, v!/’ills<^ and Bainard [ 14 ] 
found that capillaries as large as 0*3 mm diameter would gradually clog 
up under the flow of oils and would close completely if a small percentage 
of oleic or stearic acid weije added to the oil, IThis clogging was 
attributed to the gradual building up of adsorbed films on the walls of 
the capillaries, Bulkley [ 15 ] in a carefully conducted series of 
experiments with capillaries of approximately the same size as those 
used by ¥illson and Barnard [14J , was however able to show that the 
clogging was due to the fact that the oil had not been thorou^ly 
filtered. From further experiments with hi^ly filtered oils, Bulkley 
concluded that there is no rigidity in clean liqtxid at distances of 
0,5 X lO”*^ ram from a solid boundary. This, althou^ in accordance 
with the two experiments of Bastow & Bowden [ 16 , 17 ] j 'who detected no 
evidence of non-U ewtonian behaviour upto a distance of lO""^ mm, is 
in c oitradiction to Griffiths [ 18 ] , who claimed to have found evidence 
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of rigidity in hydrocarlDons of low molecular weight with gap i-etween 

10 to 2 X 10 mm. Results similar to this have also been obtained 

by Tausz and Szekely [l9] } who found that a film of lubricating 

011 under a steel ball dipped into mercury was remarkably stable, 

Ron-polar oils gave a film about 5000 A (52; 10"*^ mm) which broke only 
after tvro days, Films of polar oils maintained themselves for as long 
as sixty five days, Terzaghi [ 20 J found the viscosity of water to be 
many times the bulk value in the passage between glass plates 10 mm 
apart, Macaulay [ 21 J also found a similar increase in viscosity. 

Of great importance and interest is the work of Reeds [ 22 3 
in 1940, who squeezed fluids between two steel plates and found that 
a stable residual film was formed, Reeds* critics expended most of 
their energy on deciding that the fluids used were contaminated with 
solids, Reeds’ defence partially allays this situation at least for 
olive and castor oils , Though one must view Reeds’ work with some 
reservation, it appears very reasonable that an alteration in the 
liquid film occur ed making it substantially more viscous than in 
the bulk. 

In 1949 Hennikor [ 8 ] published an extensive 
review paper which included 174 references, largly experimental ones, 
in support of his conclusion that the surface zone of a liquid is not 
merely a monomolecular layer, but a multimolecular layer. More recently, 
but much more briefly, Adamson [ 23 ] in I960 admitted that although 
all the evidonoe of deep surface orientation is more or less circumstantial. 



10 


nevertheless deep molecular orientation does occur, prohahly as a result 
of perturhations passed on from molecule to molecule. 

•^n important class of problems v/'here the evidence of a greatly 
enhanced viscosity in thin films close to a solid surface is strongest, 
is when small proportion of a surface active compound is added to the 
fluid. In fact, some maintenanco practices reccmmend the addition of a 
small amount of used oil to new oil that may be placed in a lubrication 
system. The improved performance due to used mineral oil, as boundary 
lubricant, is probably duo to some oxidation of the oil and the formation 
of long chain acids in the oil. The advantages derived from adding small 
amounts of long chain acids to a mineral oil were recognised as early as 
1920 by Wells and Southoombe [ 24 ] , 1 little later in 1921, Hardy 
and Doubleday [ 25, 26 ] have reported tho effect of chain length of 
additives in the lubricant, They observed a definite relationship between 
tho coefficient of friction and the length of moleoulo ; longer tho 
molecule lower tho ooefficiont of friction, Beeok ot.al [27 J , who 
used a very slow running four ball machine, showed that a coherent 
oil film, separating the surfaces, was promoted by the addition of a 
small proportion of oleic acid to white oil. 

Perhaps tho most telling support for Needs observation is that 
provided by Puks [28-29 ] , Ho performed a number of exporiments with 

apparatus Very similar to that used by Needs, The resemblances to Needs* 
data are striking. It is of interest that an addition of o.l?^ of stearic 
acid to a ohromatographically separated naphthene-paraffin fraction oil 
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(]yE-20) increased the residual film thickness. It is equally interesting 
to note that lower molecular T/eight hydrocarbons such as benzene, hexane, 
cyclohexane, isooctane and decane, exhibited no residual film, Fuks 
investigated many properties of the film and it was observed that film 
thickness increased with fatty acid concentration, with surface energy 
of the solid substrate and with increasing hydrocarbon chain length 
of both solvent and fatty acid. 

The works of Needs [ 22 ] and Fuks [ 28-29] have been a matter 
of great controversy and criticism for a long time by many Trorkers, most 
recently by Hayward and Isdale [ 50 ] • It been suggested that 

asperities on the solid surface might account for the residual film. Also 
dirt and undissolved material in the fluid have been invoked as 
explanations. According to Dmuglis et,al L 51 ] (who conducted a series 
of experiments v/ith utmost care and observed residual films similar 
to Needs and Fuks), if the data of Fuks are all to be explained on 
the basis of dirt, then one must postulate dirt which changes its | 

properties with temperature, length of additive molecule, concentration | 
of additive, but such a dirt would be a most unusual naterial, | 

Derjaguin L 52-33 J , who is regarded as the leading modern j 

exponent of the view that molecular orientation in depth occurs near I 

1 

! 

solid boundaries and gives rise to significant rheological effects, 1 

j 

devised an experimental method in order to avoid any possible interference 

by dirt particles. Increase in viscosity was observed where traces of i 

( 

1 

polar or surface active substances are present (including ordinary lubricai 

j 

oils) * ! 
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llore recently Cameron and his coworkers [ 34 - 37 ] have also 
studied the influence of surface active compounds in various lubrication 
systems. The lubricants were pure paraffins, minly hexadecane, and the 
additives consisted of various long chain polar compounds. It was found 
that the oil film was markedly influenced by the additive. The scuffing 
load was found to increase with chain length of additive and solvent. 

Squeeze film studies similar to Puks [ 28-291 were also made. The 

-4 

plates came to rest at a separation of about 2 x 10 mm, with all 
fluids, even v;ith pure cetane. In addition, it was found that the 
increase in surface viscosity was greatest Tirhen additives and carriers were 
matched iii terms of chain length and shape. The results presented 
in these studies are most easily explicable by postulating the presence 
of a surface film some 10**’^ mm to 10 mm thick which has a much larger 
viscosity than the buUc of the fluid. 

Thus, although there may bo some controversy, as to the existence of 
rheological anomalies (like viscosity enhancomeiat) in the neighbourhood 
of solid boundaries, the evidence on whole _seeiE3 positive. It is -evident 
that such zones of esahanoed viscosity, giving rise to a residual film, 
vTOuld have important effects in the lubrication of practical machinery. 

Ordinary mineral lubricating oils contain various polar compounds, 
many of them surface active, in a non-polar hydrocarbon medium and 
thus belong to the class of liquids for which the evidence of the 
effect, viscosity enhancemoat, is strongest. The rheological anomalies 
could therefore have a considerable effect on the minimum film 
thickness, which in turn influences wear and fatigue failure of machine 


components 
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1.6 DEVELOPim^S 111 THEORETICAL BIIBOLO&Y 

A proper theoretical explanation of the residual film, reported 
by various \TOrkers, is clearly needed for a better understanding of 
thin film lubrication technology. An attempt in this direction is the 
work of Allen and Drauglis [ 9 ] » who developed a theory of boundary 
films based on the ordered liq.uid model, in which the wall forces, 
though not strong enough to bind more than a few molecular layer, induce 
long range order in the fluid. In this model, films on solid substrates 
are assumed to behave as liquids of special type, that is, they can flow 
as an ordinary liquid but have an internal order much like that in a 
crystal. This model fits qvialitatively for all of the data of Puks [28-29 j • 
A beginning of a quantitative theory based on molecular consideration 
is made in [ 9 ] pursued further in [58] • 

These detailed studies, although necessary for a fundamental 
understanding of physical and chemical phenomena involved in the formation 
of these films, are hardly justified for an engineering purpose. In 
such situations it may be sufficient to take into consideration the local 
degree of freedom in some reasonable average sense. For example, the 
local rotation and spin of clusters and their stretch, in addition to 
classical degree of freedom, associated with a continuum may provide an 
adequate next step. 

Motivated by these considerations, attempts have been made by 
many researchersto extend the range of applicability of continuum mechanics 
(particularly in the study of fluid behaviour of rheo logically complex fluids 
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in situation where micromotion plays an important role) using the field 
theorjj' approach to describe the macroscopic manifestation of microscopic 
events, i.e. microelement tflotions and defoimiations , 

The inadeq.uaoy of classical continuum approach to describe the 
mechanics of complex fluids has led to the development of theories of 
microcontinua in which continuous media are regarded as sets of structured 
particles possessing not only mass and velocity but also a substructure 
with which is associated a moment of inertia density and micro-deformation 
tensor. This extension of fluid mechanics req.uired a reappraisal of 
classical concepts (e.g. the symmetry of stress tensor, the absence 
of couple stresses, etc.) in order to account for local structural 
aspects and mioromotions . In fact, while many of the principles of 
classical continuum mechanics remain valid for this new class of fluids 
they had to be augmented with additional balance laws and constitutive 
relations. The presence of microscopic elements in a fluid gives rise 
not only to classical Cauchy stresses but also to couple stresses due to 
the microelement interactions. Further the orientation of the elements 
comprising the miorocontinuum had to be considered as well as the 
production of internal angular momentum due to their rotational motions. 

No attempt is made in the present thesis to give an account of 
the historical development of the theory of microcontinua, as an excellant 
review has been given by Ariman, Turk and Sylvester [59] j a reference 
may be made to that* 


The earliest formulation of a consistent general theory of fluid 
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microcontima is attributed to Eringen. [ 4 O- 4 I ] , in which the mechanics 
of fluids with deformable microelements are considered. The basis of 
Eringen’ s theory of "Simple Ificrofluids" was his work vfith Suhubi [42-45] • 

Although these papers were concerned with a nonlinear theory of microelastio 
solids, their treatment of motion, balance of moments, conservation of 
energy and entropy production is applicable to all continuous media 
where microelements are considered. 

In these papers [ 42-45] the authors develop a physical 
model in which each continuum particle is assigned a substructure, 
i.e. each material macrovolume element contains microvolume elements 

which can translate, rotate and deform independently of the motion of 

i 

the macrovolume ; however each deformation of the microvolume element 
can be expected to produce a subsequent deformation of the macrovolume 
elements. Thus a mechanism is provided in the theory to treat materials, 
whic-h are capable of supporting local stress moments and body moments, 
and in addition are influenced by the microelement spin inertia. 

The field equations (i.e. balance of momenta, conservation 
of mass and energy etc.) of Eringen* s micro continuum mechanics have 
been developed by applying Cauchy* s law of motion and the classical 
conservation statements of mass and energy in tho forms of volume 
and surface integral operators to each microelement contained in a 
material microvolixme element. Simple statistical average of these 
global balance and conservation laws have been taken over the entire 
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macrovolume element, resulting in the fundamental la^/'s of microcontinua. 

The laws of classical continuum mechanics have heen augmented with 
additional equations, which account for conservation of microinertia 
moments and balance of first stress moment which arise due to consideration 
of microstructure. Ihe key point to note in the Eringen’s theory of 
microcontinuum mechanics is the introduction of new kinematic variables, 
e.g. the gyration tensor and microinertia moment tensor, and the addition 
of the concepts of body moments, stress moments and microstress averages, 
to classical continuum mechanics. 

The micro fluid theory [ 40 ] is however too complicated and 
the underlying mathematical problem is not easily amenable to the solution 
of non-trivial problems in this field. This led Eringen [44] to 
postulate a subclass of these miorofluids, called micropolar fluids, which 
would still exhibit the effects arising from particle micromotion but 
would produce considerable mathematical ease. Under these assumptions 
deformation of the fluid microelement is ignomed, Nevertheless, microrotatic 
effects are still present and surface and body couples are permitted. In 
this theory two independent kinematic vector fields arc introduced, the 
vector field representing the translation velocities of the fluid 
particles and the vector field representing the angular (or spin) velocities 
of the particle (i.e, microrotation vector) , 

Many theories similar to Eringen* s micropolar fluid theory have 
been advanced by various research workers in recent years from different 
view point. However the resulting field equations can be reduced to 
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thoso of Bringen [44] (for a detailod comparison sgg rof. 39) • 

Physically, the theory of micropolar fluids may serve as a 
satisfactory model for description of the flow "behaviour of polymer 
fluids and fluid siispensiona, where micromotion may play an important 
role* Another class of problems where micropolar fluid theory may 
provide an answer is' when the wave lengih of an acoustical distur"bance 
becomes comparable to the mean grain-size or distance in a body. In 
such situations diffraction patterns appear which can not be reconciled 
on the classical mechanical grounds * Under these circumstances the 
molecular or granular constituents of the medium are excited individually 
and the intrinsic motion of the material constituents mast be taken 
into account. This situation could be thought of as one, representing 
the phenomena associated with thin films in proximity of a solid surface. 

1.7 mcRocouTBTumyi analysis of fluid flows 

The theory of micropolar fluids has been applied extensively 
to describe the non-Newtonian behaviour of real fluids such as colloidal 
solutions, fluid suspensions and sometimes blood, A review of all such 
published work till 1974 made by Ariman et.al, [45] j 'wl'-o listed 
over a hundred references. Since then scores of other papers have been 
published and no attempt is made to list them in a chronological manner. 

This section is restricted, basically to the discussion of those 
flow problems, which give evidence of rheological effects in the flow 
of micropolar fluids through narrow passages . 
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Increased viscosity due to the presence of microstructure in 
the fluid have been reported by Kline and Allen [ 46-48] ■ Cowin [ 49 ] 
and coworkers [ 50-51 ] » while studying the polar fluids [ 49 ] (a- fluid 

identical to Eringen^s micropolar fluid), have analyzed a special class 
of flows, called creeping flows. The most significant result of their 
v/ork is that the polar fluid exhibits a resistance to motion, greater 
than or eq.ua 1 to that of the viscous fluid in the traditional flow 
situations. Willson [52] and Peddieson and MoiTitt [53] have also 
demonstrated increased effective viscosity in boundary layers in micropolar 
fluids . 

The origin of rheological abnormalities is then to be investigated 
on the basis of the theory of micro-mechanical motions, as there is a 
distinct possibility that micropolar fluid theory, which takes into 
account the intrinsic motion of the material constituents, may provide 
answer to some of the questions concerning these thin films, at least 
in a macroscopic sense. 

Motivated by these considerations lubrication theory for 
micropolar fluids has been developed by many authors. The earliest 
work in this direction is by Allen and Kline [ 54 3 » who, using their 
own model, gave an approximate solution of the slider bearing problem . 
Similar theories were developed by Agrawal et.al, [ 55 ] » Balaram and 
Shastri [ 56 ] $ Shixkla and Isa [ 57 ] ^-h-d Prakash and Sinha [ 58 J 

using Eringen’s [ 44] model for micropolar fluids. 

Application of these theories has been made to various traditional 
lubrication problans. The slider bearing problems have been solved by 
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Datta [ 59 i t Maiti [ 60 ] , Shukla and Isa [ 57 J j Isa and 
Zaheeruddin L6l] f Agrawal [62] and Sinha and Singh l 65 J * 

Prakash and Sinha [58,64,65 ] , tehanti [ 66 ] and Isa and 
Zaheeruddin [67] have solved the problem of one -dimensional 
journal hearings. The squeezing flow in parallel plates and curved 
plates have been solved by Agrawal et.al, [ 553 , Balaram [68] , 

Prakash and Sinha [ 69-7o3 , Ramnaiah and Dubey [TI] , Zaheeruddin 

and Isa [72] and Tjyidon and Jaggi [ 75-74 J • Ihakash and Sinha [75 3 
and Sinha [ 76 j have studied the dynamically loaded infinitely long 
journal bearings. Hydrostatic and externally pressurized bearings 
have been analyzed by Khader and Vaohon [ 77 3 » ^hukla and Isa [ 78 } 
and Isa and Zaheeruddin [ 79 3 • Prakash and Christensen [ 80 ] and 
Sinha [ 8l3 have analyzed the problem of rolling contact bearings • 

The elastohydrodynamic zone of roller bearing has been analyzed by 
Prakash and Christensen [sz] . The inertia effects in squeeze 
bearings and thrust bearings have been studied by Mahanbi and 
Hamanaiah [ 85 3 • 

The earlier theories were based upon the assumption that the 
Newtonian viscosity and micropolar viscosity coefficients are 
constant across the lubricant film. Recently Tipei [ 84 ] has 
given a more generalized form of Reynolds equation which assumes 
the viscosity ooeffioioits to be function of all co-ordinates. 

However, Tipei [ 84 3 too, has solved problem of short bearings 
under the assumption of constancy of viscosity coefficients. 
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The results obtained by using micropolar fluid as lubricant 
are encouraging and illustrate the capabilities of this theory in 
explaining the phenomenon of the erihanoement of effective viscosity 
in thin films. Hotfever, the results in general are grossly 
exaggerated due to the inadmissible values of the governing parameters. 

Even after a decade this theory seems to be in a preliminary 
state. Scores of lubrication problems have yet to be solved from 
this view point. For example, little work has appeared on dynamically 
loaded bearings and porous bearings from the micro continuum view point . 
Moreover, the problems of finite dimensional bearings and bearings 
with rou^ surfaces from this view point are yet to be analyzed. 

Hence, in this thesis, an attempt is made to study these 
lubrication problems from the microcontinu-um point of view. 
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CHiPTER II 

BASIC jE!QUATIONS AlCD GENERALIZED RBYIIOLDS EQUATION 
FOR mCROPOTAS FLUIDS 

2 .1 CTTR ODUGTIOH 

The theory of microfluids introduced hy Eringen [ 40 ] deals 
with the class of fluids which respond to certain microscopic effects 
arising fi»m the presence of microstructure and are influenced by spin 
inertia. This theory was later simplified to the theory of micropolar 
fluids [ 44 j ill which the deformation of microelement is ignored. These 
fluids contain a suspension of rigid particles with individual motion 
which support stresses and body moments [ 40,44 j * This subclass is 
obtained from the general microfluid theory by imposing skew symmetric 
property of the gyration tensor, in addition to a condition of microisotropy. 

In view of the inadequacies of the classical Newtonian theory, 
lubrication theory for micropolar fluids has been developed by many authors 
[ 54»55»56,58 ] for two dimensional bearings, and by Shukla and Isa [ 57 J 
and Tipei [ 84 ] throe dimensional bearings , Apart from the analysis 

of Prakash and Siiaha [ 58] , for the two dimensional case, no paper gives 

a detailed order of magnitude study for the three dimensional problem. 

Hence in this Chapter, a detailed order of magnitude study is made and a 
generalized type of Reynolds equation is obtained for the three dimensional 


case , 
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2.2 BALMCS EQUATIONS FOR MCROPOLAR FlUIIS 

The ha, lance equations governing the heha-'/iour of micropolar 
fluids have been derived hy Bringen [ 44 ] and the results are simply- 
presented here, 

(i) Field Equations 

The field equations for micropolar fluids in vectorial form 
are [ 44 ] 

•|| + V .(PV) = o (2.1) 

(X+ 2U) vv.v _ -i (2p+)0 V ?c (Vx V) +XV»<v-vn + pF = p (2,2) 

44 ^ 4 / 4 4\ 44 ^ 4^ TT^ / \ 

(s+e+y) VV . V ~yV.}< (VxW -tX^ ^ F - 2X V + pi « p 3 — (2 .5) 


The first equation is the principle of conservation of mass 
and the others are the conservation of linear and angular momentum, P 

4 

is the mass density, V is the velocity vector, v is the microrotation 

^ 4 ’,^.. 

velocity vector, H is the themodynamic pressure, F is the hody force 
per unit mass, L is the body couple, j is the microinertia constant, X,w 
are the viscosity coefficients of the classical fluid mechanics and X,“,^,Y 
are the new viscosity coefficients for micropolar fluids. 


denotes material differentiation. 

For incompressible fluids V.Y = o, and n is replaced by the 
hydrodynamic pressure p, to be determined from the boundary conditions. 
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Hie quantities X,u,x»aj3>Y j are constants and characterize 

the particular fluid under consideration. It may he emphasized that 

the constant X links together, through these field equations, the 

velocity vector Y, and the microrotation vector v. If mass, time 

and length are denoted hy M,L and T respectively, then ^,y, and X are 

all of dimension while andy are all of dimension 

2 

and j of dimension L . 

(ii) Constitutive Equations 


Constitutive equations for the stress tensor t^^ and the 
couple stress tensor mj^^ ,are given in cartesian co-ordinates as 

*1^ - (- n* + + („. |x) 


^ ,1c ” ^ k2r 




r,r kA 


k,A A,k 


( 2 . 4 ) 

( 2 . 5 ) 


where is an alternating tensor. An index followed hy a comma 

represents partial differentiation with respect to the space variable 

(iii) Boundary Conditions 


The boundary conditions are 

Y(Xg,t)=Y^ (2.6) 

v(Xg,t) = ( 2 . 7 ) 


where Xg is a point on the solid boundary having prescribed velocity Y^ 
and prescribed microrotation velocity .The zero spin or hyperstick 
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condition, Covfin [49] and Alien a/nd Kline [ 54 ]> re pres on bo the 
extreme case in which the wall acts to completely inhibit microrotation. 
The zero spin boundary condition permits recoyery of classical results 
by passing to the limit as a suitable measure of siibstruoture size 
tends to zero. 

(iv) Three Dimensional flow 


For subsequent study, attention is restricted to three 
dimensional steady motion of an incompressible micropolar fluid under 
the following assumptions 

(1) Body forces and body ooLiples are iiegligible 

(2) all the characteristic coefficients are y-independent . 


Thus the velocity, the microrotation velocity and the pressure 
have the forms 


Y = (u(x,y,z), v(x,y,z), w(x,y,z)) 

V = (v^(x,y,z), V 2 (x,y,z), ^^J(x,y,z)) 


p * p(x,y,s) 


The balance equations (2,l) - (2.3) reduce to 
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( 2 . 11 ) 


( 2 . 12 ) 
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(2 .14) 


(2.15) 


( 2 . 16 ) 


(2.17) 


2 *5 LUBRICATION PROBIM-ASSUMPTIOI^S MD BASIC EQUATIONS 

To state the luhrication assumptions a three dimensional slider 
type of hearing is shown in Pig. 2*1 as it is the simplest load carrying 
mechanism and an easier geometry to visualize the properties of more 
complicated films, U is the velocity of upper plate along ssi^iirection 
carrying a load W, h is the film thickness | a2is the x-dimension and hgis 




Slider 
Oil film 

— ■»> X 

Bearing 


Fig. 2.1 Slider bearing configuration 
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the z ~d.imens ion of the plate. The fluid is an incompressihle micropolar 
fluid . 

The lubrication films to he analyzed are assumed to comply with 
usual assumptions i. 85j86 j made in deriving the Reynolds equation 

1. Plow is laminar ; no vortex and no turhulsnoe occur anywhere 
in the film. 

2. Body forces are neglected, e.g. the gravitational or magnetic 
field is assumed to he negligible, 

3. The film is thin enough as compared to the length and span 
of the hearing. This permits us to ignore the curvature of 
the fluid film and the rotational velocities may he replaced 
hy translational velocities. 

4 • No slip on the hearing surfaces • 

Apart from the above mentioned conventional lubrication 
assumptions, three more assumptions must he emphasized. Later on, in 
Chapters VI, VII and VIII, the first two will he violated, 

1, Bearing surfaces are smooth, i.e. there is no surface 
roughness * 

2 , Bearing surfaces are non-porous • 

3, The hearing surfaces are rigid, i.e. they are not deformed 
under pressure caused hy fluid film. 

Equations (2 .12) - (2 ,17) can he simplified with the help of 
above mentioned assumptions and following non-dimensional scheme. 
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X ^ J Z i- 


X = - , y 


a ^ 


h 




^2^ 


U 


TJ 


w 


V .h 
1 


y 1- 


P = 




Ua2(P‘+ ^ X) 


6 — 6 — 2i 6 — 

^ 1 r% ^ 0 ^ 




2 "bo 3 ar 


^ = f ^ " ^ 2 yix ^ ~ 


The equatioiis (2,12) ~ (2,17) ill diinensioiiless_ form are 
9x 3y 9Z^ 3y 8 z 


Ip t, r - 3u-aa,^-3tt 

= Re {6-|U + v — ;;;■ + 6'2^“7r 

3x ■*• 3x 3y 3Z 


„ „2- ^2- 0.2- - 3v^ Sv, 

, 2 a_T^S^ 52 lv 1_ i 

^ 35= 3y2 ^3i2 sS 3® 


1 T, rx — 9"’?’ “ ix “ t 

ii = {6 u — - + Y -— +6 2^ — - } 

1 By Sx 3y ' Bz 


r X 2 9 D 2 9 r X i::. ±\ 

^ 3x2 gj2 2 j 52 1 -r. .r- 


2 ^^^2 


3x ?y 


. . Ee {5 ,u -^ + T ^ + 6 /r ^ } 

^3 si 9x 3y 32 


2 - 2 - 2 - 

3v, 9v, o^v-, 

{6 ? i + i + 5, } + 




1 -=2 ’ ay2 ^ gg S' 2(l-M^) - “ ^ 


ax-^ 3y 


. Sj -^} 

3y 9? 


3 Vq 


o — Vq « Re {6,u 

(i-T) ^ 


3X 


+ Y + ^2^ 

9y 


3Vq 


( 2 , 18 ) 


(2.19) 


(2 , 20 ) 


(2 , 21 ) 


} ( 2 , 22 ) 
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2 2 — 2 
O 3 V o 3 Vo ’^o 

+ + {2 _^J + ... ” i .fs -2a _ . -2H_, 

^ ^”2 2 g -2 2 ( 1 - K ^) ^ ^ ~- 


3y‘ 


3z 


8X 


- ^2 




V2 = Se’ {6^u — - + 


9x 


{6 




1 .-S ■ ,J2 • “2 ..-J" ^ 


dV^ 

3y 

3y 


9v, 

C 

9Z 


+ -} ( 2 , 25 ) 


au 




} 


9Z 2 2(l-K^) ^ 9x 9y 

3v -7 SV'T 3v ' 


V- . Re' <6,s -:::i . v 

(1-N^) 5 ^ 9X 




9y 9 Z 


} 

(2.24) 


where 


Re 


2 P hU 
2V-+X 


Re' 


P.iUh 

2 

AvS- 


(2.25) 


The parameters N and i characterize the micropolarity of the 

fluid and distinguish a micropolar fluid from a Newtonian fluid. The 

parameter H, called coupling number, couples equations ( 2 , 2 ) and ( 2 , 5 ) , 

For SLnd these equations become uncoupled and (2,2) reduces 

to the usual Nccrier - Stokes equation, Theimiodjmamio restrictions 

2 1 1 

require that o R < 2 [ B4»87 J • 


The parameter £ has the dimensions of length and can be thought 
of as a fluid property depending upon the size of the fluid molecule. 

In this sense the parameter L( - h/£) characterizes interaction between 
the geometTcy and the fluid. The case L “ (in which the substiruoture 
size is negligible as compared to the bearing oleamnce) gives Newtonian 


results 
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It may be noted that Bringen [ 44 ] , Shnkla and Isa [ 57 ] 
and Tipei [ 84 ] and other authors have used the parameter NL (where 
L = h/i,) as one parameter, in discussing the micropolar fluids, hut 
this product characterizes a combination of effects, which are considered 
separately in this thesis. These parameters are discussed in more 
details in next Chapter. 

Ee is the modified Reynolds number where u has been replaced 
by V + ■^ X and this number is always smaller than the classical Reynolds 

X 

number. For lubrication problem Re '^■^ 1 , generally of the order of 10 '« 

! 

Also it is reasonable to assume that Re << 1, since 3, is the square 

2 

of the length typical of the microstructure and oon^iarable to £ , If 
a requirement that miororotation velocity identically equal to the 
local fluid vortioity should be a possible solution of the balance 
equation, is imposed, it can be shown that 3 = £ [ 47 ] • This requirement 
seems to be a natural one and has the desirable effect of reducing the 


number of constitutive coefficients. 

-5 

Further, 6^ and fig are of order 10 , 6 ^ is of order unity. 

Using these order of magnitude study and the above requirements, eqns . 
(2,19) to (2.24) reduce to 


lii + 2l!l^ 

3y 2 




( 2 . 26 ) 


PI = 


O 


(2.27) 
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a 

21^ fi 

1 

-^ = 0 

CM 

CO 

ay 

■*“5 


2- 

^ji 


X . 

3w 

3^2 

2(l-N^) 


ay ( 1 - 1 ^) 

. 2- 




-Ji __ 



ay 2 

l-N^ 

K V2 

- 0 

2_ 

B Vj 


an , 


3y2 


ay 

(L-1^) 


(2 . 28 ) 


(2.29) 


(2.50) 


(2 .51) 


Prom eq,uation (2,27), P = p(x, 2 i), arid when (2,30) is solved with 
zero boundary conditions [ 49>54] » ^^2 

(i) Velocity Distribution 


To solve for velocity components and microrotation velocity 
components, we revert equations (2,26) ~ (2 . 51 ) to corresponding 
dimensional form . 


The equations are 

_ „2 9v, 3p 

(P+ ix) ^ + X — =0 

ay ay ax 


-1 ^2 av- 

(„+ iy') .y i » = 0 

IP+ ^^2 9y az 


2 ' 
a^vi 


^1 o-,/ 9''®' 

Y 2Xv_ +x -r— = o 

2 lay 

ay 


(2.52) 

(2.55) 

(2.54) 
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9 ^V3 

Y ^ - 2xv^ 

sy^ 5 



o 


(2.55) 


The four differential eq,uations (2,52 to 2,55) form two 
systems of simultaneous equations. One consists of (2,52) and (2,55) 
and the other of (2,55) s,nd (2,54). The solution of first system will 
yield expressions for u and The second system will yield expressions 

for w and , it is seen that the second system "becomes identical to 
the first system if we replace "by -v^. Hence it is necessary to 
solve only the first system for u and and the expressions for w 
and can "be written in a similar way. 

The "boundary conditions for velocity and m icrorotation 
velocity are 

®= w = ii^ 2 » Vp ~ °f "^3 = y = 0 (2,56) 

^ " ^21’ ^ “ ^22* '^1 “ "^5 “ y = ii (2.57) 

The solutions are 
2 

~ p ^Tx 2 “ ~m ^"^21 "^51 my} + ■^■^ 2 ( 2 . 58 ) 

^5 “ iiT'^Tx ^ "*11^ "^21 “^51 

’ = p( ^ 2 + h2=^) - < Ajg sinh i»y + cosh ny}+ A42 

(2.40) 

Vi = - { ~( y + Aj^) + ^22 ^ "^52 (2-41) 
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where i “ > j = 1»2 | are constants and given by 


[ (lJ„ . - IJ .) sinh mh ih sinh mh + “(l - cosh mh) } ] /A 

J 


A, . 

A = 

2j 2y 

A,. = i[ {n„.-u .) d , r °°Bb . ■■h) + (cosh * -X) 

53 y 23 13^ 2u 8x. 2 " 

^ 3 


(2 .42) 
(2 45) 


+ (h ~ ~ sinh mh) } ] /A^ 


^43 


2iNr^ 

Ut . + — A, . 

l3 m 53 


2M^ 


5 “y 


{ sinh ^ (oosh mh « l)) 


(2.44) 

(2.45) 

(2.46) 


in which 


X. 




:^ = x, X2“Z»ni = Y 


(2.47) 


(ii) Flow Flux 

The volume flow rates and along x and z axes, are 


given by 


Sc- 1 


(2 48) 


and 


q_ = / w dy 
2 -'o 


(2.49) 
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or 


2 ^ li 2V \1 9X 


(2.50) 


and 


~ 2 ^^12 ^^^22 


) _ h^f(H,Ah) 3P 
y 3 z 


(2.51) 


where 


2 

f (N, £,h) - 2 + , 2 “ 2h 2£ 

h 


(2.52) 


It should be noted that for cornrenience in Chapters VI and VII the floi» 
fluxes are defined as follows 


^x 2^^11 ^21"^ y 3x 


== hfjj +p ) _ £ ( ^^1 ^ ) _ii 

O\^0^ '^OO ' %, 


"■z 2^ 21 "22 
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(2.55) 

(2.54) 


where 


\ h^, , .2 H £ h ^ ,, Kh 

f(l!f,£,h} = -^2+ h £ - -g — ooth ^ 


(2.55) 


(iii) Generalized Eeynolds Equation 

Integrating the eq.uation of continuity (2 .ll) across 
the film and using (2,50) and (2,51) (or eq.n 8«2.53 2,54) 

together with (2,56) and (2,57) the generalized Reynolds equation 
for mioropolar fluid in three dimensions is obtained as 


' “ '>■ 


9x 


+ p + 7 _V 

21 9x 22 9z oh 


(2.56) 
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where and are the normal velocity components of the surfaces 
at y = h and y = o, respectively. 

Similar equation was obtained by Shukla and Isa. [ 57 ] 

and Tipei [ 84- J under the assumption of constancy of characteristic 
coefficients across the film. However the parameters H and ~ 
are different from the parameters used by [ 57»84 ] • 

For infinitely long bearings, along z-axis, the two 
dimensional Heynolds equation is given by [ 58j • 


(iw) Frictional Torques 


The shear stresses aloi^g x-direction and z-direotion, 
from equation (2 #4)} 

- («+ + XV, 


X 


(2.57) 


and 


\ -if - 


(2.5S) 


and on the surfaces y = o 8.nd y = h these assume the form 


( T ) , 

x^Ojh 


- h ^ ^^2l’"^ll^ 

"2 gx hg(H,£ ,h} 


(2.59) 


(t ) 

^ z'o^h 


- h 9P '^^^^22"^12^ 
■^2 gz hg(N,Jl,h) 


( 2 , 60 ) 


where 


g(H, 2.,h) 


h ^ ^ 2£ 


( 2 , 61 ) 
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The frictional dre^s are obtained, by integrating the 
shearing stresses over the bearing surfaces. 

It should be noted that for convenience in Chapters VI 
and VII the shear stresses are defined as follows 


\ = 2 li — 

^ x'ojh 2 3x g(Ni2. ,h) 


, . - h ap 

*^\''o,h 2 32 g(N,£,h} 


( 2 . 62 ) 

( 2 . 63 ) 


g(N,J!,,h) = h ~ 2WS, tanh ^ 


( 2 . 64 ) 


where 
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GMITER III 

TWO 11‘D THEEE DEISt:STO?liL FOF-C'rCTIG SOFgEZE PILIffi 
3.1 IMTEODUCTIOH 

In recent years isany workers have analysed the problems of 
micropolar fluid lubricated non - cyclic squeeze films [ 55j68-71,83 J 
and cyclic squeeze films [ 75-76 ] . The problem of porous non-cyclic 
squeeze films have also been solved using micropolar fluid as lubricant 
[ 72-74 } • The prominent feature of using micropolar fluid as 
lubricant was an increase in load carrying capacity and squeeze time. 

Though some of the geometries, such as circular plates, 
considered were finite dimensional, but due to non-variation of 
pressure along circumferential direction, the problems solved were 
essentially two dimensional. It is thus imperative to study the 
effects of micropolarity on three dimensional squeeze films problans. 

In this Chapter, an application of micropolar fluid theory 
is made to study some hitherto untouched two dimensional in^ortant 
problems and. to some three dimensional configurations. 

The problem considered is that of the so called ’ non-oyclio 
squeeze films* . The squeeze velocity V is imparted by a load W, which 
will be considered constant in magnitude and direction and symmetrically 
placed, with respect to the boundaries of the system, Tarious geometries 
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are considered , namely f spherical bearings, slider and rectangular 
plates, conical bearings and some miscellaneous configurations. The 
approach is similar to that of Pinkus and Sternlioht [85] • I^etailed 
numerical results are obtained for the sjiierical bearings, whereas 
for the three dimensional case only the dimensionless response time 
is obtained. 


5.2 SPHERICAL BEAHINGS 


The film geometry for squeeze film in spherical bearing is 
shown in Pig, 3.1(a) , Here we shall consider the case of hemispherical 
seat . 

Prom the geometry, given by Pig,5.l(a), the film thickness 
is given by 


h =» r - R - a^cos e (3.l) 

which can be rewritten as 

h >= o(l -e cose) (5*2) 


where 


<= = r -K| ^ - ^3- 

Por a sphere the amount of fluid passing through a conical element 
of stir face is given by 


(3.5) 


h 

q a / (SttR sine)u dy 

^ r\ 


(3.4) 


where u is the velocity of the fluid, Por zero tangential velocities 
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the expression for u, equation (2 *38), reduces to 


u 


oosh my - In h dp 
"yEdS^f” m s ihh mh 2iiE d9 


y “ "T ^ 
m 


X {3ihh my - 

sinh mh 


(3.5) 


in which x is replaced hy Ee * 

Substitution of u, equation (3.5) equation (3.4), yields 


q^ = 2Tr E sin0{. f(N,il,h)} 


(5.6) 


where m,E and Z are given hy eqn« (2,18) and f(E,£,h) is defined 
hy eqn, ( 2 * 52 ) . 

The displacement of the fluid at any point in the seat is [ 88 ] , 

q^ = TTYE^sin^e (5.7) 

where V is the relative normal velocity of the surfaces, responsible 
for squeeze film action and is given hy 
de 


7 => a 


dt 


(5.8) 


Equating (5.6) and (5.?) 

2 

dp ^ uTR sin 9 

2h^f(E/,h) 

The pressure distribution with the boundary conditions 


(3.9) 


p( + 2 ) “ o 


( 5 . 10 ) 


becomes 
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ir 

.2 7 


p(0) - < 

0 ,h) 


ir 
1 - 2 


( 5 . 11 ) 


since the prescure distrilmtion is synimetric for the regions - ^ 
and. 0 ^0 <_^ . 


ilon-i'-iraencionalizing, using 
H = ™ = (l_ e cos 0 ) , L = ” , P = 


me 


2o 


( 5 . 12 ) 


the non-dimensional pressure is 

ir/P 


P(e) = / 


siri0 d0 


0 h5p(it,l,h) 


( 5 . 15 ) 


where 


P(1T,L,H) = 


1_ 1 _J__ 

(LH)2 “ 


coth 


KIH 


( 5 . 14 ) 


The load capacity for a hemisphere is 


2 * Q . -r^VlR^ sin^0 dl 

= Ptt R j p sin 6 cos 0 de = — J 

n d. ^ 


O h^f(ir,il,h) 


( 5 . 15 ) 


The non-dimensional load is 


W 


2o^ 


^ sin^e.de 


0 H^P(¥,L,H) 


( 5 . 16 ) 


The biaie of approach, t, ftom an initial film thickness 
h^ to a final film thickness \(l\^ terms of cojcres ponding 

values of eccentricity ratio can he evaluated using ( 5 . 15 ) and (5.8). 
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Thus 

■'r/2 5 

j ' ^ sin'^ 9 d9 de 

0 0 ■ h5f(lT,il,h) 

The noru-dimens ional time of approach is 

rp _ t _ j sin^ 6 d6 de 

(.liiSi.) S ° H^F(¥,L,H) 

2o^¥ 

3.5. SLIDERS AM) RECTMGTJIAR PLATES 

The rectangular plate poses a three dimensional prohlam, 
and we must return to the hasio differential equation of a viscous 
flow of mioropolar fluid for a solution* 

The differential equation governing the pressure distrilmtion 
in three-dimensions, equation (2*56), for h = constant, is 

(5.19) 

3x 3Z h^f(N,^,h) 

Eiquation (5*19) is Poisson's equation, the solution of 
which is assumed of the form 

CO 

p(x,z) = \ (A cos Z(2i) ' (3.20) 

n=l,5,5 ®'2 

where 2(z) is the function of z alone and is a oonstanb to he 

determined is the x-dimension of plate. Substituting (5.20) into 
(5.19) and using the relation 


t = 


frUoR"^ j 
e 


2W 


(3.17) 


( 3 . 18 ) 
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“ (n»l) 

1 ^ ^ 
n«l,3»5 


1 mrx 
— cos ““ 
n 


we have 


I (ps.f A cos Z( 2 ) + 

n=l,3,5 ti ^2 


I A n^x d Z(z^ 

i.. cos 

u«l,3,5 


^2 dz^ 


= i 


UV 


n„5 


I . , (-1) 


(»=ll 


h^f(ll,^,h.) «,=1»5»5 


1 niTx 

— cos 

tl &n 


which is satisfied if 


A. 


Asm. (1^)2 A^zw-f^ Tcvh?- 


(n^ll 

2 


dz^ 


A particular solution is given "by 

nrrz „ , nro L’^XL- 

Z(^) = sinh — + 0^ ocsh ^ 

Putting ( 3 .J( 1 ^ into ( 3 - 20 ) toundaiy oonaitions 

\ 

p »s 0 Bjh 2 s= + 2 

where hg z-ddnension of the plate and hy the requirement 

symmetry in z-direotion, the pressure distrihution is 

33^1 


nrrz 

a 


HU ” (-l") ^ 

-I3.5 A 


- 1 


cosh 


2a, 


(3.21) 

(3 .22) 

(3.25) 

ll 

(3.24) 

(5.25) 

of 

] (5 *26) 
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The total load is given by 


¥ = + 


16 V a^Y 




« a iiTTbp bp 

^ ^ “5“ “2ap “ Tlf 

n=l,5>5 T'’' 2 2n^ 


) 


(3.27) 


dh 

Using time of approach from an initial film thickness 


h^ to a final film thickness h^(h^> h^) is given by 


h. 


dh 


16 p a^ “ b a„ nv b ""1 

—T I ( — 4 - “1^ ■fcanli / T7 

n=l,5,5 2n^ xPk ^ h^f(U,Ji,h) 


(3.28) 


Kon-dimensionlizing, using 

h 


h 

h. 


H . ¥ , H . 




T 


( 


16 V &2 


J ( !h tanh 

n.1,3,5' 3^4 - 


nvb , 
2a„ 


(3.29) 


the non-dimensional time of approach is 

-T = /^ --JS (5.30) 

H^P(N,L,H) 

3.4 ELLIPTICAL PLATSS 

For an elliptic plate 



(5.31) 
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the pressure distribution with the help of governing differential 
equation (3 .19) can he written as 


2 2 

p(z,z) = ( -^ + — - l), C is a constant 

^ 'P,‘~ 

1 1 


(3.32) 


where the pressure boundary condition is p = o on the boundary of 
the ellipse given by equation (5 .31). 

Substituting equation (3.331) into (3.19)j '''’e ge't the 


pressure distribution 


p(x,z) 


2 2 


2 2 


2(A^+B^)h^f(N,ii,h) "^l \ 


(3.33) 


The load capacity and time of approach from an initial film 
thickness h^ to a final thickness h^ are given by 


W = + 


\h 




(5 .54) 


^ uA^B^ h 

1 ^ ^ f dh 

(A^+B^) h^ h^f(N/,h) 


(3.55) 


Circular plates are special case of elliptical plates and the bearing 
characteristics can be obtained directly f3X>m above expression by 
putting A^=B^=r, Prakash and S inha [ 69 ] • 


5,5 MISCELLMEOUS COM' ICBRATIOKS 


By the methods similar to those described above, results for 
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other shapes, of falling plates, can be extracted. Some of these 
are given by 


(i) Circular rings of inner radius r^ and outer radius 

2 2n2 

u2-*i-L ^ 


\7 


IT yy 


r „4 r ^^2 " ^1^ 1 

8h5f(N,£,h) 


(3.56) 


putting 


_ M 


Y = - 


dt 


t = 


Try ^ A A ) 

8W" ^2"-"l' 

log(#) 

1 


2 2x2 h 




dh 


h h^f (f,£ ,h) 


(3.37) 


(ii) Circular sector of angle B j^ and bounded by radii r^ and rg 


¥ 




Tr^h^f(N,Jl,h) 


4 4 
2 1 


I [ 

^•5>5 4n^2--al-)2 

®1 


2Tr 


1 


1 


O+SIL 2+SI- 

r2 -r^ 

2niT 2nTT 

6-1 B- 

rg 1 - 1 


( 5 . 38 ) 


and 


8 


t ~ 


I 

n«l,3,5 


4 4 
2 1 


,nrr 


C- tQ 

r„ B 


2 "1 - "^1 '^1 


2 + 


Sol 


. 2,^ nTT .2 

4n {2^j — ) 


2njT 


2mT 

r ®1 
^1 


2tt 


n6^{4-(f^)^} 


dh 


/ 

h h^f (N,a ,h) 


X 


(3.39) 
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Letting o in equations (5.38) and (3 *39) the characteristics for 
full sector can he obtained, 

(did) circle of radius r . ^ , a portion of which has been cut out 
by another full circle of radius r 2 (fig, 3.1(b)) 


Tf 


MY 


r 4 , 2 2 An -1 

L (2r!r-4r r -r •) cos 


8h^f(h,2',h) 


■i^-"r2 -^2 


2r, 


4* 


+ ^(4+ I r|) (4x^-4) 1/2 j 


(3.40) 


and 


* " + ^ 2^ + *2 4 '!'' 


( 4 r^- 4 )V 2 ]/"_^ 


h h^f(lT,«-,h) 


Letting rg o , the case of fall circle is obtained, Prakash 
and Sinha [ 69 ] * 

(dv) Full Conical Seats 

In case of complete cone (Pig,5 . 1 ( 0 )) , 


¥ 


Trp 

•-t- 


8h^f(N,i!<,h) 


(3.41) 


(3.42) 


where ds the semi-vertical angle and is the radius of base 
and 

Try e:; 


h 


t = + 




dh 


SWsin^a^ h^ h^f(N,il,h) 


( 343 ) 
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Fig 


3 1(c) Full conical seat configuration 
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(■'^) l^runoated Conical Seats 


For truncated cone (Fig, 5»l(d)), whose snaller radius 
is Rg and larger radius is R^, 


W 


TTUY 


8h^f(lT/,h) sin'^u^ 


p Op 

r r4 r4 1 

L R2-*2“ r -* 

log(^) 


and 


t = 


■'P .,4,4 dh 




8Wsiji\ (hs ll„ hh{ii,4,h) 

^2 

The dimensionless response time in each of the above 
oases is governed by integral 


/ 


dH 


H- h5p(R,L,H) 


This integral is therefore integrated numerioally for 
rectangular plates and represented graphically. This integral 
is similar to the one obtained by Rl^akash and Siriha [ 69 ] • 


3' ,6 RESULTS MR DISCUSSION 

(i) Dimensionless Parameters 

The mioropolar fluid is characterized by three physical 
constants, U ,x and Y as compared to a Newtonian fluid which only 


( 3 . 44 ) 


(3.45) 


(3.46) 



J I 



Ftg 3 1(d) Truncated 


conical seat configuration 
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viscosity, y =Uy + [44] • Thus the expression for II hecomes 

p-2 . __X 

2(y^+x) 

Thermodynamic restrictions require that , X ^o [ 40j44] » 
it is thus evident that o ^ . 


(3.48) 


The second parameter L, defined as the ratio of the clearance 
width to the characteristic material length ^ , characterizes interaction 
between the geometry and the fluid and may he considered to he a property 
depending upon the size of lubricant molecule, say, the chain length 
of the molecule. 


It is expected that the effects of micropolarity will he of 
importance either when the characteristic material length is large, 
corresponding to a large size of lubricant molecule or when the 
clearance width is small. The second case is of importance in 
lubrication theory where o is usually small and hence comparable with ^ . 
Recalling that L = o/l and keeping c fixed an increase in I means a 
decrease in L, so lower the value of L, larger the chain length 
of molecule and higher should be the effects due to micro polarity. 

In fact it is seen that as L ->■ o, the expressions for velocity and 
other flow characteristics reduce to their equivalents in the Rewtonian 
theory, with d everywhere replaced by (u + ^X) , Thus indicating an 
increased effective viscosity. This observation supports the experimental 
evidence of increased effective viscosity t 8,22,28—36 ^ in the presence 
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of additives having large chain length and as -well as in the case of 
thin films. Since a small value of L corresponds also to a small 
clearance width provided £ is kept as constant , A large value of L, 
on the other hand, signifies a large clearance v/idth, i.e. a case of 
hydrodynamic lubrication or a small material characteristic length 
and the effects of the micropolarity are not likely to he significant. 
It is seen that as L tends to <» , £ would tend to zero, and the 
individuality of the microelement is lost and once again the Uewtonian 
results are obtained, 

(ii) Bearing Charaoterist ics 

Out of the various cases that have been treated in this 
Chapter, the numerical results are obtained in detail for the case 
of spherical bearings (Pig, 5*2 to 3«7) whereas for the three- 
dimensional rectangular plates only the response time is plotted 
(Pigs , 3 •B and 3 •9) • 

Non-dimensional pressure distribution curves are shovirn 
in Pig, 3*2 (N = 0,7 and various values of L) and Pig, 3 . 3(1 = 0,0 
and various values of N) , Since the pressure cuive is symmetric 

/ TTv 

about 0 = 0 , only half of the pressure curve is shown (.0 <_ 6 . 

L “(Pig. 3,2) and N o (Pig, 3*5) show the Newtonian pressure 
curves. As the value of L is decreased, which corresponds to a small 
clearance or a large size of additive molecule, the pressure begins 
to rise (Pig, 3 •2) and for L = 0 , 0 , the naxifflum pressure 

is obtained , As the coupling number increases 
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(Fig. 3*3) the additives play a significant role. The pressure in all 
cases is higher than the Newtonian pressure. 

Fig. 3,4 and Fig, 3 *5 show the graphs for dimensionless load 
capacity as a function of micropolar parameters N and L» The load 
also is seen to increase as the value of N increases and L decr^ses. 

The Newtonian curves are straight lines as the load remains unaffected 
by the variation in N and L in this case. 

Figs. 3.6 to 3*9 the curves for response time for various 
values of N and L. Response time is actually the most important 
squeeze film characteristic. This gives us the time that will elapse 
for an oil film to be reduced to some minimum value, a value which 
may have been established as the thinnest permissible film for a 
particular bearing. The time required to reduce the film thickness 
can then be compared with the length of the time the load is actually 
applied. 

Fig, 3*6 shows the time-height relationship for a fixed 
value of N(N = 0 . 7 ) and various values of L, For large values of H, 
i,e# for largo film thickness, the effect of additive on response 
time is not significant . But as the surfaces further approach each 
other, the value of H decreases and a marked discrepancy is observed 
between the Newtonian and micropolar curves. This discrepancy is conspicuous 
for H < 0.3 and for low values of L, A similar result was obtained by 
Needs t 22 ] in his experiments on squeeze films. This discrepancy 
indicates that fcr low values of L and H, the time of approach is 
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increased, whioh. in turn indicates the eriharcement of effective 
viscosity, which is a confirmed experimental fact '.dien additives 
are used in lubricants. Similar interpretations can be made 
regarding fig. 5»7» 

Pigs. 3»8 and 3*9 depict the response time curves for 
various values of L and N for rectangular plate squeeze film 
(three-dimensional problem) . The qualitative effects of micropolarity 
on response time is same as that in two-dimensional problems which 
can be easily ascertained from these graphs. 

It can thus be concluded that the theoretical effects of 
additives on three-dimensional lubrication problem are almost 
identical to those in two dimaisions, at least qualitatively. 



180 



Fig. 3'2 Non dimensional pressure distribution for 
various values of L at N = 0-7 and £-0-7 
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|.3‘4 Non - di mcnsior-al load vs. L for various value 
at e=0-7 


of N 









62 



T - — > 

Fig. 37 Non-dimensional height vs. non-dimensional 
time of approach for various values of N 

at L-2-0 


UJ ; 



time of approach for various valuer^ of L 
at N = 07 (Rcctar.gular plate) 




Fig. 3-9 Non-dimensional height vs non-dimensional 
time of approach for various values of N 
at L = 2*0 (Rectangular plate) 
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CHAPTER rV 

DYHAI, lie ALLY LOADED SHOHT JOTJEIAL BEAEINGS 

4.1 IHTRODUCTIOH 

Dynamically loaded bearings in which the load alternates 
or rotates has been subject of much interest, conjecture and research. 
Dynamic situations where both wedge and sq,ueeze films are active, arise 
in the bearings of almost all reciprocating machines, especially in 
steam engines, automobile and aircraft piston oigines and various 
machine tools. Even the human joints can be approximated by 
dynamically loaded bearings* 

The load-«arrying phoiomaaon in such bearing arises due to 
the fact that a viscous lubricant cannot be instantaneously squeezed 
from between the two approaching surfaces. Because of the resistance 
of the lubricant to extrusion, pressure is built up and the load 
is supported by the oil film. When the load is relieved or reversed, 
the oil film thickness can recover its thickness before the next 
cycle, if the bearing has been designed to permit this build up. 

This phenomenon is found in reciprocating machines where tl® 
bearings are subjected to fluctuating, impulsive or dynamic loads 
and therefore present problems not met in the normal operation of 
steady load bearings* In general, lubricating oils in the bearing 
subjected to reciprocating loads are hi^ly vulnerable to contamination, 
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Toy iDoth dirt and metal pa-rticled. The lubricant may then hehave as a 
fluid suspension. Moreover, the use of additives is a common practice 
not only for the enhancement of viscosity and a consequent increase in 
the load carrying capacity, hut also to increase the working range of 
temperature. These additives, through chemical reactions, with the 
metal surfaces form solid surface films and thus protect the hearing 
stirfaces from damage and wear. 

Little work has heen reported [75»76] on the luhrioation 
of the dynamically loaded hearings from this view point. It is thus 
felt that the application of the concept of mioropolar fluid lubrication 
to hearings operating under dynamic conditions would lead to a better 
understanding of hearing operation and could yield results in conformity 
more to practice than those obtained by the Newtonian theory. The 

i 

application of the present work may also he found in the lubrication 
of human joints, which can he approximated by dynamically loaded 
hearings. These joints are lubricated by synovial fluid, which has 
long chain hyaluronic acid molecules and can thus be thou^t of as a 
mioropolar fluid. 

In earlier works, by Prakash and Sinha [ 75 ] and Sinha [ 7^ 1 » 
the problems solved were restricted to infinitely long journal bearings. 
However, for practice, the short bearings are better approximation to 
finite journal bearings [ 89 ] ♦ Hence, in this Chapter, we solve the 
problan of dynamically loaded short journal bearings lubricated with 
micropolar fluid. 
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¥e start with the general case of dynanic loading in which the 
load and speed of rotation vary hoth in magnitude and in direction with 
time, Various hearing characteristics for an infinitely short journal 
hearings are obtained in the form of integrals, Hie equation for 
deteinining the journal locus is rather complicated and therefore only 
a particular case is considered in detail, name ly^ pure squeeze films 
with no journal rotation. The case of a fluctuating (sinusoidal) load 
without rotation is considered. Under this loading condition, the 
journal may he thought to act like a piston in hydraulic damping 
mechanism. Physically this implies that, as the load is applied 
in one direction, the fluid is squeezed out on one side and sucked in 
on the other. Before the piston comes in contact with the hearing, 
the load reverses its direction and the oil on the other side begins 
to take the load, Tims, assuming a continuous fluid film to exist 
from o to 2Tr the hearing response to sinusoidal loading is studied 
and the mioropolar effects are elaborated graphically, 

4.2 THE PROBmt nil) MAIYSIS 

The problem considered is that of dynamically loaded journal 

hearings (Pig,4»l) in which the load is a function of time. The journal 

is in combined tangential motion U and normal motion V^, where U and 

V are hoth functions of time, 
o 

Equation (2 ,56) , governing the pressure distribution, for 
the type of flow mentioned above, assumes the form 
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3 




{ ^"f(N,£,h) —} 
az p ^ ^ ^ 9z 


3P 


2 -12i + li -il + Y 

2 3x 2 3x o 


(4.1) 


where N and £ are given hy ecLn. (2,18) and f(N,^'jh) is defined "by 
eq.n, (2 ,52) , 

The first term in the r.h,s, of eqn, ( 4 .I) represents the action of 
the journal rotating with a velocity U over a wedge shaped fluid film 
given hy h(x), the second term implies a variation of tangential 
velocity along the hearing surface and the last term is the expression 
for velocity of the shaft centre and is responsible for squeeze film 
action. 


Assuming the hearing to he fixed, the shaft centre will 

have instantaneous radial and tangential velocities. Referring 

d G d (4> +^) 

to the vertical line ip = o, those will he c tt-* and e — ^ 


dt 


dt 


Thus U and are given hy 


IJ = Rm +0 sin 8 ~ cg 


cos 6 


( 4 . 2 ) 


and 

Y = O' cos 0 + c e sin e (4.5) 

o dt dt 

Substituting for TJ and Y^ in equation ( 4 .I) and considering 
that, L85J , x=R6,|«2 and ~ cose« 2 and 
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sine- 

dt dt do 


( 4 . 4 ) 


one obtains, the Reynolds eq.nation for the general case of dynamic 
loading for mioropolar fluids, 

i- 9 rh^ f(n.J^,h) 1 + J_ { h"f(R,£.h) , 


R‘ 


,2 80 


89 


8 z 


3z 


r A, o,> n ^ 4 h-, , „ d£ Q 

i -2 ^} + 2o ^ cos e 


( 4 . 5 ) 


For short journal hearings, eq.n. ( 4 . 5 ) reduces to 


dz 2h^f(lT,A,h) 


^ { (to. - 2 ^ + 2o,^ cos 0} (4.6) 


/ (iif; 1 

■where co is the angular velocity of the journal and = "j:^) is the 
frequency of rotation of applied load. 

Integrating eqn, ( 4 . 6 ) twice, using the boundary 
conditions 

L 


p( i ~) - 


( 4 . 7 ) 


the pressure distribution is obtained as 


p(e,z) = 


4 h^f(ir,fi',h) 


{ (a>-2 £0 t- 2 'l^) ~ + 2o cos 9 } (z^-7^) 


dt 


(4.8) 
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The load component, , normal to the line of centers is 


given hy 


^o 

T 

= ¥ sin^= 2 / / p(9,z) sinSR de dz 

2 o o 


(4.9) 


Sulostituting the value of p and integrating, we get 


W = ¥ sin4> 

TT 


2^ . 2 . .a 

Of ^ ^ d$ xr sin 6. d9 

-(a,- 2(0, - 2 — )i “ 


24 


o h^f(lT,£,h) 


(4 . 10 ) 


The load component, W^, along the line of centers, 


is given hy 

A 

2 2it 

W = W cos (j> = -2j f p(e,z) cose d6 dz 
° ■'o o 


(4.11) 


which after simplification becomes 

cy L^R - 2^ ^^2 

•r^r ^ d£ f 003 ^ ay 

Q TO ’" aT ^ J ^ 


12 dt o 


( 4 . 12 ) 


Equations ( 4 .10) and ( 4 .12) give a relation between 
the applied load and the motion of shaft center. In these equations 

W, (0 and oan be functions of time. This makes the integration 

’ L 

of these equations more complicated, bat, wten accomplished this 
would, through the relation between e and 4>, srield the cyclic locus 
of the shaft center, as well as the instantaneous resultant 
hydrodynamic force . 
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(i) Cyclic Sq,ueeZ3 Films in Journal Bearings 

In. this section, the stndy is restricted to the particular 
case of dynamic loading, in which the journal does not rotate about 
its center hut moves under a variable unidicectional load. The load 
.vector is fixed in the space, and the shaft center follows a path 
determined by the nature of imposed load. 


Here , since 


0 ) = 0 )^ = o 


we have from eq.ns . (4»10) and (4 •IS) that 

5 

vecL R , . ZT . 2 . 
W, =Wsini.= — 


n 

2 


12 dt ^ h^f(N,)l,h) 


and 


3 

PCL-^R 2^ 2- _ 

W = W cos ^ 

° 12 dt J h5f(H,£,h) 


(4.15) 


(4.14) 


(4.15) 


Dividing, we get 


de 

d^ 


tan<j) 


fl(^) 

'^(iJ 


where 

and 


sin^e d9 
o h^f(H,Jl,h) 

cos^Q d6 
o h^f(N,A,h) 


(4.16) 


(4.17) 


(4.18) 
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Integrating (4 •16) , we get 


s in (j) = A 


exp[/ 


d£ 


] 


where A is the consta.nt of integration and 






(4 .19) 


(4.20) 


Equation (4»19) describes the locus of the shaft centre 
under a unidirectional load, which depends upon the cost ant of 
integration A, determined by the initial position of the .iournal. 
Equation (4 •19) is independent of load which shows that the journal 
centre once located on a given path will always remain there regardless 
of the magnitude and nature of the applied load. The effect of 
load will only be upon the amplitude and speed on that fixed path. 

The substitution of (f) from eqn,(4.19) into either of the equations 
(4 •IA) and ( 4 « 15 ) would yield the relation between the load and 
eccentricity, 

(ii) Bearing Eesnonse to Sinusoidal Loading 

We now consider the case of shaft travelling along the 
vertical centre line(4'= o), with an oscillatory load, 

W(t) = sin (4 •Si) 

where W, is the amplitude of applied load and a is the frequency 

^ ir 


of oscillation 
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Equation (4 •15) becomes 


sinupt 


ycL^ _ 2T^ 2 . 

o de r COS 0 de 


/ 


o h^f(ii,£,h) 


( 4 . 22 ) 


■which after simplification becomes 

■fTcsy 


de 

dt 


12 W1 sinm t ^ 
_b _ 1 


(4.25) 


UEoLq 2 


Using the non-dimensional scheme 


H = li = 1 +e cos e , L = 7 , T = u)’ t = 
c ^ A ’ p 


60 


720W 

p 1 

o L D 


(4.24) 


and defining the Sonmerfeld number, S, for an infinitely short 
bearing as 


_ un - ['^\2 

“ P '^0 ^ 

0 


(4.25) 


eqn, (4 •23) reduces to 


dT 


sinT 


PS f,(2p) 


(4 .26) 


where, D = 2R, and 


2P 2 


5 N / co s 0 d9 
^2(2^) -/ P(U,L,H) 


(4.27) 


and 


F(H,L,H) 


12 


^ “2111 2 


(4.28) 
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Equations (4.26) and (4.2?) are integrated numerically 

to yield the time variation of the eccentricity ratio over the 

region of T from o to 2'rT , The resulting response curves for 

various values of ITjL and S are shown in Figs. 4*2 to 4.10. The 
L 2 

value of (^) is taken to he o.o5. 

4-4 RBSULTS MD SISCIBSIOh 

One of the most important hearing characteristics is the 
response time, i.e. the time that will elapse for an oil film 
to he reduced to some minimum value, say a value which may have 
been established as the thinnest permissible film for the safe 
performance of a particular hearing. 

Figs. 4,2 to 4,4 illustrate the journal centre locus 

for various combinations of the micropolar parameters 31 and L and 

the Sommerfeld number S, The value of is kept fixed 

L 2 

((■g^) = o.o5). For smaller values of N and larger values of L 

the effects are not significant • . But for higher values of 31 
and smaller values of L (where the micropolarity is significant) 
the effects are quite significant. This shows that the micropolarity 
causes the fluid to offer more resistance to journal motion, thereby 
allowing for smaller eccentricities (higher fihn thicloiess) for a 
constant load. 

The velocity of approach, de/dT, is plotted against time, 

T, in Figs. 4. 5 to 4,7, over a half cycle. It can easily be seen 
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that the velocity of approach decreases for mioropolar fluid for 
those values of T which are near to Tt/2, hat for those values of 
T which are close toT=oorI'=Tr, the trend is almost same* 

The velocity of approach is highest for a hewtonian fluid hat 
for mioropolar fluid this velocity decreases. Thus indicating 
an increased effective viscosity. 

The maximam eccentricity ratio, e , is a measure of minimum 

' max’ 

permissible film thickness, and thus its study is essential. The 
maximum eccentricity ratio is plotted against L in Fig, 4«8, against 
S in Pigs. 4 *9 Fig« 4»10» In. these figures it is seen that for 

larger values of H or smaller values of L, the maximum eccentricity 
ratio decreases. For the Newtonian case it is highest. It implies 
that the same amount of load is now being carried with a smaller 
eccentricity (or at a larger film thickness). This means that the 
system can he used to sustain much higher loads, or, conversely, 
the same load can he used for a longer period of time without any 
aal func t ioning • 



Eccentricity ratio 





4.3 Path of the journal centre o\)<er one cycle of sine 
loading for various values of coupling number 






Eccentricity ratio 








Fig 4. 7 Velocity of journal centre over one half cycle of sine 
loading for various values of coupling number 




Maximum eccentricity ratio - 



L > 


Fig 4,8 Maximum eccentricity ratio vs. length ratio for various 
values of coupling number 







Fig, 4.9 Maximum eccentricity ratio Vs Sommerfeld number for 
various values of length ratio 
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CMPTER 7 

LIJBRIGATION OP ROLLER BSiRTJC-S IIT COilBHIED 
ROLLIIIGjSLIBIRG AED NORMAL MOTION 

5 ,1 INTROBNCTION 

Thousands of papers on film, luhrication have appeared 
since 1886 in the wake of Oshome Reynolds classical v/ork on 
lubrication theory. In general, papers on film lubrication have 
extended the application of the fundamentals, as set forth by 
Reynolds, to a variety of bearing shapes with steady and transient 
films of variable density, viscosity and temperature. 

In recent years a lot of attention has been directed towards 
the study of rolling element bearings. Considerable understanding 
of the bearing characteristics has been obtained through the studies 
of elastohydro dynamic lubrication, and for all practical cases 
theoretical solutions have been published. 

Most of earlier work in lubrication have dealt with steady 
lubricating films in which there existed e'ither an effective 
tangential surface velocity or a maximum film velocity due to an 
externally pressurized source. Unfortunately there are few instances 
where the film are completely steady. Usually there are transient 
or periodic forces or displacemaits imposed on the films. Changes 
in surface velocity result in tangential accelerations which can 
be important during starting and stopping. However, relative normal 
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surface motions are of far greater importance. The film must maintain 
sufficiently precise spacing in the presence of disturhance, Various 
authors have examined theoretically the hearing hehaviour under normal, 
sliding or rotating motion soparafcely. Sasaki et.al [ 90 ] were 
the first who studied the dynamic hehaviour of roller hearing. However 
the boundary condition applied at the trailing end violates the 
continuity of flux, Howson et. al, [ 91 later on, obtained a 
detailed theoretical solution of lightly loaded cylindrical hearings 
subjected to a combined rolling, sliding and normal motion applying 
a more realistic boundary condition at the film mipture point, A 
satisfactory agreemaxt was found between theory and experiment [ 92 ] . 

No attempt has so far been made to anaHyze this problem [ 91] 
from the microcontinuum view point • Thus it is felt that the present 
analysis may provide a better insight into the lubrication of such 
bearings when analyzed from microcontinuum point of view. Howson’ s 
et,al, results [ 91 J obtained as a special case of the present 

study. It will also be shown that this theory may yield better 
fit to the experimental curves obtained by Marldio and Howson [ 92 ] , 

5 *2 THE PROBLEM MH THEORETICAL iNALISIS 

The problem considered is that of the lubrication of lightly 
loaded cylinders in combined rolling, sliding and normal motion, 
with micropolar fluids. The fluid is assumed to be isoviscous. 

The geometry of the system and the coordinate system are shown in 
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Fig,5«l.The two cylinders of radii and rolling with, peripheral 
velocities and respectively, and approaching with a noimal 

velocity V , are separated by a lubricant film of thickness h. The 
conventional representation of the system is a plane and a’ cylinder 
of radius R. Moreover, the circular cylinder can be approximated 
by a parabolic cylinder with the film thickness h given by [ 95 ] 

2 

... , .X 


(5.1) 


where 


1 1_ 
R " Rj_ 


+ 


2 ^ 


where h^ is the minimam separation of the oil film between the 
surfaces , 


(5.2) 


Assuming the cylinders to be infinitely long in the axial 
direction, the governing equation for a two dimensional flow field, 
from eqn. (2.56), is 

~ {h^f(R,Jl,h) - S + V} (5.5) 

where R and £ are given in (2,18) and f(lir,£,h) by eqn. (2,52). 
Integrating eqn. (5.5)j using the condition 

^2 

^ = o, St X = x^ or h = h^ = h^ + ^ ( 5 . 4 ) 


the pressure gradient is obtained as 
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Condition (5.IO) is the Swift-Stieher honndary condition applicable 
when I q.l < 1 [ 91 j . 


The other boundary condition depends on the oil flow 
through the T/edge. If the oil supply is enough and the weight of 
the oil is n^leoted, the pressure build up will begin at the start 
of the wedge, "where X = , This is the inlet boundary condition 

for the pressure, i.e* 


P = o at X = -00 (5»ll) 

Integrating eqn. (5*9) fcr the pressure distribution, we 


have 


X (2q4X4-Xj (XJC ) dX 
P(X) = / 2 S 


h^p(h,l,h) 


(5.12) 


The point X can be determined using P(X ) = 0, i.e. 
m m ^ 

X (2q4X+X ) (XJC ) dX 

i ? == 0 (5.15) 

H^P(ir,L,H) 


Solving eqn. (5*13) we can find the values of X^ for different 
values of the parametorsh,! and q. 


The load capacity per unit length is 


iv 

W . f p dx ~ - f x-gdx 

00 


(5.14) 


or, the non-dimensional load capacity is 
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2%“ 


u(Ui+U2) E»^2Bh^ 


= - / 




^ dX ® 


(5.15) 


The shear stresses acting on the plane (y = o) and cylinder 
(y = h), in regions where oil is filled, are given hy ein» ( 2 , 59 ) j 

hdE- 

"''Ojh "" ~ 2 dx hgClI^I jhJ 


( 5 . 16 ) 


where g(K,JZ.,h) is defined hy eq.n, (2 ,61). In the oavltated region 

the mean shear stress is [94] 

_h y(u^-h,) 

•^o = 


(5.17) 


The surface tractions are 


dx 


h,2 - ij-2to+ h^gdl.l.h) 


( R -Ir) 


or, in non-dimensional form 



21^?. 9 
0 


y (U^+U2)>^o 



— CO 


+ J 


dX 

J hg(u,l,h 5 ’ 


H / 

•m ^ 


m 


X 


m 


dX 


(5.19) 


where 


g(x,l,h) 


1 


~ tanh 
hti 


Nffl 

2 


2^ 

J - Uj_+U2 


(5.20) 
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These surface tractions are put in the form 

C^] (5.21) 

f2 = -[B^„jC^l (5.22) 


(Sl+X+X^) (X^J d X 

^ -CO 2H^p(lf,L,H) 

X ^ 

f ^ dx r dX 

m 

The coefficient of friction on the plane is 


(5.23) 


(5.24) 


h P^ 

(r)"f- '5.25) 

The numerical solution consists of first solTing the 
eq.uation (5.13) to get the values of X^ for different parameters j 
N, L, and q.. The equation is solved using Regula falsi method 
within an appropriate accuracy where the infinite limit has heen 
replaced by a finite but large qxxantity, (i.e. is replaced by -lO) . 
A comparison is made between the solutions for the steady state case 
(q = o) with the noiu-steady case (q o) . Por this purpose some 
ratios for different bearing characteristics are defined. 


(i) Location of the Pilm Bupture or Cavitation Point (X^) 
Let 

X 

X = m 

" K\.c 


(5.26) 
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where represents the general point in which q. can take any 
accepted values and is the point of cavitation for sliding 

or rolling case. For the same minimum film thickness, X is a 
function of micropolar parameters IT and L, and the velocity 
parameter q,, 

(id) Location of Line of Maximum Pressure 

It has been already pointed out that if X^ represents the 
point of cavitation then the expisssion for determining point of 
maximum pressure X , say, can he found as 

X = - 2q. - X ( 5 . 27 ) 

P 

and thus the nature of this line would he same. Defining 




(X 


m \=o 
P 


(5.28) 


where q, = 0 represents the steady state case, 
function of 1T,L and q. 


X 




will he a 


(iii) The Maximum Pressure 


Let 



P 

m 



q==o 


(5.29) 
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7/here is the value of maximum pressure under pure rolling 

or s liding case for the same instantaneous value of the minimum 
film thickness , i. e* » This is also seen to he a function 
of Kj L and q.. 


( i'v) Instantaneous Load Parrying Capacity and Surface Tractions 




and 


- ^ 

^ ^0 



(Vq=o 


(5.30) 


(5.51) 


(5.52) 


where the term ( )q^_Q X'epresaats the steady state case, then it 

F and are dependent on 
^1 ^2 

N, L and q. 


can he seen that the functions , 


The numerical solutions have been given for the point of 
cavitation, the instantaneous load carrying capacity and surface 
tractions. To have a better insight of the present analysis the 
pressure distribution, load carrying capacity and coefficient of 
friction are presented in non-dimensional foims for different values 
of N, L and q and compared with the Newtonian case (N = o or L-^"). 
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5.3 R33ULTS AlJJj I‘ISCI55I0If 


(i) 


Dimensionless Parameters t 


-^p8.rt from the micropolar paramebers F a,nd L, one more 
important parameter is encoun'*'ered in the present study, i,e. the 
velocity ratio parameter, cl, defined hy 


q = 


Y 





( 5 . 55 ) 


which determines the qualitative behaviour of a roller bearing when 
subjected to a combined rolling, sliding and normal motion. Positive 
values of q indicate a aittiation when the surfaces are moving apart 
normally and negative values signify an approaching case. Numerical 
solutions have been presented for ~o.9 < q < o.9V since, it has been 
pointed out by Dowson ^ $2 ] j that the cavitation occurs for [ q| <1. 


(iii) Bearing Characteristics s 


Non-dimensional pressure distribution curves are shown in 

Nig. 5.2 for different values of velocity ratio paramsfcer q and the 

micropolar parameter N for a fixed L, Mioropolar curves are shown by 
dotted lines, ^rom these curves it is seen that the approaching normal 
velocity (q negative) results in a stronger pressure generation while 
the separation (q positive) leads to a weakKC pressure. The 
mioropolarity ofthe fluid increases the pressure at corresponding 
points as compared to the Newtonian fluid (N = o.o). 
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Son-dimensional load capacity W is plotted against the 
length parE,meter L for different coupling numbers S and velocity ratio 
parameter q. in Sig, 3 * 3 * It is clearly seen that as L increases 
the load capacity decreases and its value is seen to ax^proach the 
Newtonian one for different values of q. For the fixed value of 
L the load capacity is higher for larger values of N and it is 
least at N = 0.0 (which corresponds to the Newtonian case [ 91 ] 
for the corresponding value of q) , A similar conclusion can he 
dmwn regarding the parameter q. For negative values of q the 
load capacity is higher and for positive values of q it is smaller. 
When the normal motion disappears (q = o) , we get the case of pure 
rolling and sliding. 

The variation of the load parameter ratio ?d.th the modulus 
of the dimensionless velocity ratio q is shown in lig. 5.4, for 
different values of the micropolar parameter N. This figure depicts 
that tho effects of micropolarity aro higher for higher values of |q| • 
The parameter in case of normal approach, increases as q and 

N increase and decreases as these decrease (Fig, 5»4)» The figure 
also shows that the load capacity is higher when the surfaces approach 
each other and it increases as the rate of approach increases or the 
micropolar parameter N increases. In case of separation tho phenomenon 
is reversed* 

The frictional drag parameters Fj^ and Fj^ , 


as a function of 
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N, L and q have heen drafra. in Fig-, 5 » 5 » interes'L'ing feature 

of these curves is the graph for F,, ‘vhen the surfaces are either 

Eg 

in normal appj:oeoh or normal separation. It is seen tiiat the 

parameter Fg^ first decreases and then increases when the sircfaces 

are moving apart for all values of micropolar parameters. But v/hen 

the surfaces are in normal approach the function Fg increases 

2 

slowly. The use of micropolar fluid is seen to reduce this factor. 

The qualitative "behaviour of the parameter Fg is the same as that 
of load parameter. 

The coefficient of friction (E/h^)p^ is plotted against 

the length parameter L for different values of the coupling numher 

II in Fig. 5,6 (q = -0,5), Fig, 5,7 (q = 0.0) and Fig, 5,R (q = 0,25) , 

The value of ^ , m ■ is taken to "be 0,5 in all these three 

1 2 

cases. It is seen that the coefficient of friction reduces to the value 
predicted "by the classical lu"brication theory for all three limiting 
oases (L o, Le-co, andll-vo). It is seen tha.t as the coupling 
number h increases the departure of tho coefficient of friction 
from the classical theory becomes more conspicuous. The most 
interesting aspect of these cizrves is the fact that tho coefficient 
of friction decreases rather sharply as L increases from 0 to 
approximately 10, The minimam value of is a function of L and 
I'T which can be seen directly from the graphs. This observation is 
true for the case of noimal approach and sliding or rolling. It can 
be also seen that as q decreases from +0,25 (Fig, 5»8) to -0.5 (Fig, 5, 6) 
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for fixed value of the coupling number H, the mininum values of 
coefficient of friction are obtained for smaller values of the 
length parameter L. Fig. 5*8 is the graph of coefficient of friction 

in case of separation of surfaces, peculiar behaviour of coefficient 
of friction is seen for this case, i-s L increases the coefficient of 
friction first increases and attains a maximam value for a fixed 
coupling number and then decreasesrapidly and reaches a minimum 
values (a value which is always lesser than Newtonian one) and then 
increases again slowly, The behaviour is more significant for 
higher value of N, 

Fig, 5,9 is the graph for point of cavitation vs, the 
velocity parameter q.. It is seen that as q. decreases the point of 
cavitation increases. The point of cavitation for pure rolling or 
sliding case is obtained at q. = 0 , 0 . In case of micropolar fluids 
the cavitation process is accelerated, i.e, it starts a little 


earlier 



4 0H\ ' 


q=0.50 


N = 0 6 


N=0.0 '-\ 




> 
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Fig. 5 2 Dimensionless pressure distribution 
of the dimensionless velocity rciuo 
different coupling number N 
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Fig 5 7 Coefficient of friction (R/h^)|j.f vs L for different 
values of N ( q = 0 0 ) 
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Fig. 5.8 Coefficient of friction (R/Fq) pt vs. L for different 
values of N(q = + 0.25) 
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CHAPTER YI 

LUBRICATION THEORY FOR ROUGH SURFACES AIID 
ITS APPLICATION TO A JOURNAL BEARUTG 

6.1 INTRODUCTION 

The analysis of surface roughness effects on luhrioation has 
gained considerahle attention and has "become a subject of intensive 
study during the last two decades. This is due to their importance 
in the practical field, ^ its physical significance and also to the 
great theoretical interest inherent in these problems. 

Most theoretical studies on "bearing lubrication appearing 
in the literature make use of the assumption that the bearing surfaces 
are perfectly smooth. It has long been recognized that owing to 
machining limitations, this assumption is rather unrealistic, 
pairticularly in bearings, where the film thickness is smll. 

In recent years a major way of analyzing the problem has 
been the adoption of statistical methods. Two basic approaches 
have been developed namely the deterministic and stochastic approach. 

The starting point in each case has been the Reynolds epilation. In 
the deteirministic approach the surface rou^iness was accounted for 
by postulating a sine (cosine) wave or a series of sine (cosine) 
waves for the film thickness [ 95 J • More recently Tseng and Saibel 
[96,97] have used stochastic concepts to study the surface roughness effects 
on slider bearing [ 96 ] and on short journal bearing [ 97 ] • The 
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analysis was concerned v/ith the problem of two dimensional slider 
bearing, with one dimensional roughness in the direction transverse 
to the sliding direction. There are, however, experimental evidences 
[ 98 ] which often exhibit a one-dimensional type of roughness 
running in the direction of sliding, generally referred to as the 
longitudinal roughness . Keeping in line with this evidence 
Christensen and Tender [ 99-100 ] and Christensen [ lOl] presented a 
stochastic model for hydrodynamic lubrication of rough surfaces, 
Christensen [ lOl J considered two types of one-dimensional roughness, 
the longitudinal rou^ness with its striation parallel to the sliding 
direction, and the transverse roughness with its striation perpendicular 
to the sliding direction, and obtained modified Reynolds eq.uations by 
making two heuristic assumptions about the pressure gradient and flow 
rate. It was concluded that the longitudinal roughness results in a 
slight decrease in load carrying capacity and an increase in frictional 
force thereby causing a significant increase in coefficient of 
friction. The effect of transverse roughness is, however, to improve 
the bearing characteristics. Later Elrod [ 102 ] in^roved the 
Ghjcistensen’ s modified form of Reynolds equation, 

' Scores of papers have appeared in the wake of the earlier 

mentioned works on various aspects of the lubrication of rou^ 
surfaces. All of the works appeared to date show that the surface 
roughness may have significant effect when the minimum film thickness 
is conparable to roughness mean height , The effect is much pronoimced 
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for lower values of film thickness. Inherent in these works was a 
fundamental assumption that all the luhrioants hehaved as if they 
Were Newtonian, 

Non-Newtonian "behaviour is almost invaria"bly observed in 
lu"brication practice and they may "be of important consequence when 
the surfaces are rough, Christensen’s models [ 101 ] do not account 
for non-Newtonian situations or where the micropolarity of the fluid 
is important, A question, therefore, naturally arises as to how 
the "bearing performance would "be modified in such situations. 

It is the aim of this Chapter to analyze tho surface roughness 
effects from micr ooontinuum view point. For simplicity it is assumed 
that the moving surface is smooth and the stationary surface is rough. 
Three cases of roughness are considered and the theory is applied 
to the pro"blem of journal "bearing under half Sommerfeld "boundary 
conditions . 

The present study may be of important consequence where the 
combined effects of surface rotiglmess and additives are desired, 

6.2 MODIFIED REYNOUS EQUATION FOR ROIHH STOlFACgS 

The Reynolds equation for mioropolar fluids for smooth 
surfaces is of the form (eqn, 2,56) 

- ir ^ 



Ill 


where 


m ^ fjiij .h) aP 

2 y 3x 


q. 


z 


JP 

y az 


( 6 . 2 ) 

(6.3) 


obtained from eqns , (2.5 5) and (2,54) by substituting , 

^21 ” '^12 ~ ^22 “ °* ^ defined by eqn. (2,18) 

and f(li,£,h) by eqn. (2,55), 


The geometry of the lubricant film can be considered to 
be constituted of two parts. The first part denotes the nominal (smooth) 
part of the geometry and can be considered to be a function of space 
and time co-ordinates. The second part of the film geometry, a 
randomly varying quantity, arises due to the surface roughness and 
is measured from the nominal level. Mathematically, 

h = h^(x,z,t) + hg(x,z,c) (64) 

5, a random variable, determines a definite roughness arrangement, 
from a large number of similar- but not identical - arrangements 
having the same statistical properties. 

Taking the expected value, equation (6.I) assumes the 

from 

-f; ^ (6.5) 


where 
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) = / ( ) f(h ) dh (6.6) 

o S 

— CO 

and is the probahility density distribution of the random 

variable h^. Using (6.6) it is readily seen [ 101 ] , that 

aE(h) ^ 

3t at 

where h^ is the smooth part of the film thickness. 

To evaluate the average of the fluxes appearing in eq.uation 
( 6 . 5 )* subject to a particular, specific model of the roughness, i.e. 
either along or perpendicular to the sliding direction, two assumptions, 
similar to those used by Christensen [ 101 j , are made . 

(i) The pressure gradient in the roughness direction is assumed 
to be a stochastic variable with zero (or negligible) 
variance 

(ii) The flux perpendicular to the roughness direction is assumed 
to be a stocbastio variable with zero (or negligible) 
variance . 

A justification of these assumptions for riioropolar fluid 
can be forvj'arded in a way similar to that of Christensen [ 101 ] • 

(i) Longitudinal. One-Dimensional Roughness 

The roughness is assumed to be running parallel to the 
sliding direction in the form of long narrow ridges and valleys. 

Thus the film geometry assumes the form 


h = hj^(x,z,t) + hg(z,g) 


( 6 . 8 ) 
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Taking the average of all terms in eq.uation (6 #2) adopting 
the first assumption the mean flux in the x^direotion is given hy 

S(q,^ « 1^3(11) - ^ ,li)) (6.9) 

This is sOj "because j sp/sx is a varia'ble uitli zero variance 
and f(l'T,£,h) and Bp/s x can "be considered (approximately) to "be 
stochastically independent q,uantities» However, to get the mean 
flux in z— direction, we divide equation (6.5) "by f (l^,^,h) . Taking 
the expected value of both sides of resulting eq.uation and using 
the second assumption 

E(n ) = » i 1 (6 .10) 

Substituting the expressions for eq.uation (6 .l) 

the modified Reynolds equation becomes 

Ij- <E(f(H,«,h)) ) + -fj < 

.^{,,2^ + 2-5321.} (6.11) 

fii) Transverse. One-d)imensional Rou^^ness 

In this model, the rou^mess striation is assumed to be 
running perpendicular to the sliding direction in the form of long 
narrow ridges and furrows . The film thickness therefore assumes 


the following form 
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h = h^(x,z,t) + lig(x,g) (6.12) 

In .a similar -wav^ adopted earlier, the flaxes, along and 
perpendicular to the sliding direction, are 




u ^^f(¥;7,hy 

2 1 N 


1 3S(p) 

t 3X 


(E(f(lT,£,h)"^)) 


(6 . 13 ) 


eCi^) = -^^^E(f(ll,£,h)) 
and the modified Reynolds eq.uation is 


(6 .14) 



(iii) Rnifoimi. Isotropic Roughness 

In this model, the rotighness is assumed to he uniformly 
distributed over the hearing surface with no preferred position 
or direction on the surface. The film thickness assumes the 
form 

h = h^(x,z,t) + hg(x,z,5) (6 .16) 

Adopting the same hypothesis as in [99 ] mathematical reasoning 
given hy Tender [ 103 ] , the fluxes and modified equations are 
given as follows . 
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3(V = I B(h) - E(f(lT,£,h) 

E(q^) = - i E(f(N,Ji,li) 


(6.17) 

( 6 . 18 ) 


~ {E{f(i;,<;,h)) { E(fC::,^,h)) -fM) 


{u + 2 

2 dx at 


(6.19) 


(i'^) distribution of Roughness Heights 

The roughness distribution function which is generally 
used to evaluate the several expected values is [ 101 ] 

f(h ) = . (b^-h^)5 -b < h < b (6.20) 

s ^2b' ® ® 

= o elsewhere 

where b is the half total range of random film thickness variable 
and if a is the standard deviation, b = + 3® . The polanaomial 

distribution function (6 .20) is an approximation to the Gaussian 
distribution. The reason for using such a polynomial distribution 
is that the Gaussian distribution always implies a finite probability 
of having asperity of very large sizes violating the conditions 
for hydrodynamic operation. 

6.3 BEARING GEAMCTEEIISTICS 

Once the modified Reynolds type equation is obtained, for 
a particular case, the evaluation of the bearing characteristics 
can be carried out in a way similar to that of the smooth bearing. 
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(i) Oil Flow 


When the roughness is longitudinal the mean fluzc along 
and across the direction of sliding is given by eq.uation (6.9) and 
(6.10), respectively. For transverse roiighness these values are 
governed by (6.15) and (6.I4). Equations (6 .17) and (6.16) account 
for these characteristics for the isotropic case. Hie total flow 
over the edge of the bearing can be found from the a-bove mentioned 
equations by integrating in a usual way, 

(ii) Load Capacity 

The mean load carrying capacity is obtained by integrating 
the mean pressure over the bearing surface, 

E(f) = / B(p) dA (6.21) 

A 

(iii) Frictional Drag 


The expression for shear stress on the moving surface, 
from eqn, (2,6 2), is 

^ h 8P , uU 


( 6 . 22 ) 


where g(ir,jj,,h) is defined by eqn, (2,64). 

For longitudinal roughness the mean shear stress is given 

E(0 - I 


(6 .25) 
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-j_ 3 3 : 

T(¥,“a,iiy^ 


(6.24) 


The mean frictional dx&g is 

bC?) = / S(t) dA (6.25) 

A 

(iv) Coefficient of Friction 


The mean coefficient of friction is defined hy 



( 6 . 26 ) 


6,4 MOhlFTTHT) rshtqids equatioit applied to a rough journal bbarius 

The modified Reynolds equation derived above will now 
be applied to the lubrication of one-dimens ion?-l rough journal 
bearing under half Sommerfeld boundary conditions. 

The configuration of the system is shown in Rig. 6.l(<t) . 
Here the bearing is assumed to be rough and journal to be smooth. 
The two cases of roughness^ i.e. longitudinal and transverssj are 


treated separately 
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(i) Longitudinal , One-*Uimens ioual Rou^hnes s 

Since x = equation (6*ll) aseumea the form 

iSM 


de 


(6.27) 


The smooth part of the film thickness is 


h^ = c (1 + e c C6 O) 


(6 * 28 ) 


and 

E(h) = h^ (6.29) 

The houndary conditions for pressure are 

e(p) = o, at 6 =0 and tt (6'*50) 


The mean pressure distribution is given by 


0 ^ 

E(p(e;) = 


where 



IT h de 
f a 

/ d 6 

E(f(N,£,h)) 


(6 .Jl) 


(6.52) 


The mean load components pe^: unit length, along and perpendicular 
to the line of centers, are obtained by integrating the mean 
pressure around the bearing surface from e=oto 0 = it. Thus the 
mean load component normal to the line of centers per unit length 
is given by 
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IT 

E(W^) = e(¥) Sin<> =/ i:(p(e)) sin0R d0 
2 ° 


(6.35) 


Substituting the value of E(p( 0)) from eqn. (6.5l) and integrating 
we get 


S(W^) 

2 


= ^ yUR^ / 


(h -h ) cos 9 d6 
n n 

0 


■B(f(iT,T;h}y 


(6.34) 


Similarly, the mean load component along the lino of centers 
is given 


E(W ) = E(W) cos (>= ■r pUR^/ 


(h -h X . 

ir ‘ n n ) sxn 0 de 




The resultant mean load capacity is given by 

1 

E(W) = {(e(W^))^ + (E(W 


(6 .55) 


(6 .36) 


The expression for mean shear stress assumes the form 
e(t) « + mjE( — - — r-) 


g(N,^h)' 


The mean frictional drag, e(f) , per unit length is 


given by 


E(F) - /' ^ ^ IS 40 +» J “ ^0 


2Tr 

J 

o 


(6.57) 


or, using eqn. (6 .31) » 
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h (h -h ) de 
e(F) ~E(f(N,Jl,h)') ' 




The coefficient of friction is obtained according to 


equation (6.26). 


(ii) Transverse. One~Dimensional Roughness 


For this roughness the Reynolds equation is 




i B TO ^ } (6 40) 

^ ^ "(-fra) 


In a way similar to the case of longitudinal roughness, the mean 
value of pressure distribution, load carrying oapaoity, frictioml 
drag and coefficient of friction are 


e(p(0)) « 


^ ^ ®^f (R,£',h7^ “ ^n_®*^fT¥rrrhT^ ^ 


n^ 'llR,5.,h) 


(6.41) 




-) } cos 0 d0 (6 . 42 ) 


:(w^) .mj 

O 


-)} sin 0 de 


(6 .43) 


•- ° ®^f (N/ ,hP 








f(R,4;h)^ 


(6.44) 
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where 


h 


n-, 


/ E(' 


h 


-) de 




(6.i5) 


The coefficient of friction is obtained as usnsl. 


(iii) Isotropic Rotighness 

For one dimensional case the isotropic roughness will 
he eiuiwalent to longitudinal roughness. 


6.5 ITON-BIMBHSIOiTAL FORMS 

To obtain the results mmerically the following non-dmensional 


scheme is adopted, 
h 


h 


h 


h 


H 


s 


h = ^ = 1 +ecos 8 , ~ , 


B 


!.£. 5(,) , HWT = 2U|Lo" . 2i(rlo 


S' ■ UUR pUR" pTIR 

( 6 . 46 ) 

The non-dimensional average values of bearing characteristics 
for longitudinal and transverse roughness are given as foi:iows. 


Longitudinal, one-dimensional roughness 


0 (H-H ) d0 

e(R(Q)) - / s (h,l,h,b) ' 

0 X 


( 6 . 47 ) 


ir cosede 

'(6 •48) 

E(W^) = B(¥) sin 4> = j 
_ 0 

“s^TwP) 

TT 

sin8 d 8 

(6 .49) 

e(W ) = E(¥) cos d = / “■ 

° n 
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, IT H(H-H ) de -n 

-(S') 2 j S '(1I,L,H b) 2 / C-(i'i,L,H,B) da 

0 1 0 

(6.50) 

1 

E(f/)' = (irw^^+ S(W^)^)^ 

(6.51) 

2 


_ Kn 

f(E(yp) iM 

(6.52) 


where 


H de 


E = 
0 


o 1 


s. (K,L,H,b') 


TT 

/ 


ae 


s^(n,l,h,b5 


(6.55) 


S^(H,L,H,B) = 


* 7 64^1 f (H« ) 


g(e,l,h,b) = / 


B 


(B^-E^)d H 


(6.54) 

(6.55) 


52B ^ ^ ^ 


Transverse, one-dimensional roughness 


E(P(0)) 

e 

= { 8 ^( 11 , L,H,B) - S 2 (E,L,H,B) } de 

(6 . 56 ) 

i(w T 

TT 

2 

7T 

= {S^(E,L,H,B) - S^CEjLjHjB)} cos e de 

(6.57) 

E(W^) 

Tf 

= / { S^(E,L,H,B) - S^CNjLjHjB)} sine de 

(6 .58) 

E^py 

IT s^(e,l,h,b) 

2 (j ^ S2(E,L,H,B) “ 



2ff S^(E,L,H,B) 

+ 2 J {5S4(E,L,H,B) + G(1,I.,H,B) -5 -^tlT^/BT 

> de ( 6 . 59 ) 
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where 


/ s ae 

o y 


r"'’’ 


/ de 


( 6 , 60 ) 


S (N,L,H,B) = 

32B' _B 




H 


[ (H+Hg)5 + 12 ^ coth ] 

( 6 . 61 ) 

To determine the rate of change of hearing diaraot eristics 
for a rough hearing as compared to the smooth theory the following 
hearing parameters are defined to study the micropolar effects. 

(iii) Mean Load Ratio Parameter 


The mean load ratio parameter is defined hy 

E(W - 




( 6 . 62 ) 


■(^(WTIb. 


B=o 


(iv) Mean Frictional Brag Ratio Parameter 


The mean frictional drag ratio parameter is defined hy 


R 


m'-jsm 


B=io 


R 


( 6 . 63 ) 


Ie(p)) 


Bs=o 


( v) Mean Coefficient of Friction Ratio Parameter 


The mean coefficient of friction ratio parameter is defined by 



126 



6.6 HESULa?S MB DISCUSSIOl 

(i) Dimensionless Parameters 

Three dimensionless parameters are of important consequence 
in context with, the present study, namely, N, L and B, defined 
hy equations (2.18) and ( 6 , 46 ). 

The parameters IT and L have already Been discussed in 
Chapter III, 

The parameter B( = '^ is a characteristic of hearing 

c 

geometry and arises due to the surface rou^ness. It is obvious 
that the hearing characteristics would he prominent when h hecomes 
oomparahle to c, detailed discussion and limitation of this 
parameter are given in [ 99t 100,101] • Here it suffices to say 

that o ^ B < 1, 

(ii) Limiting Gases 

The present study can he considered as a generalization 
of the works of Christensen [ 101 J and Prakash and Sinha [ 58] 5 
because their results can he obtained from the present analysis hy 
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considering certain limiting cases. 

Christensen' s results are obtained in the limit L “ 
or IT -V o, as both of these cases correspond to the Newtonian results , 
The results of Prakash and Sinha are obtained in the limit B 4- o, 
since B ->■ o signifies a smooth surface. 

(iii) Bearing Characteristics 

The analysis of surface roughness effects from a microoontimum 
view point, yields bearing oharact eristics which are functions, not 
only of the roughness parameter 3, but also of the mioropolar 
parameters N and L. 

The load ratio paraneter , friction ratio parameter 
and the coefficient of friction ratio parameter C^^ hare been plotted 
as a function of the roughness parameter B for Tarious values of the 
coupling number !!(£= 0.5, L = 30) ir. Pigs.6-2to 6 .4 •It is seen 
that, akin to the Newtonian lubricants, the transverse roughness 
leads to a significant increase in load carrying capacity Pig. 6.2. 
and the frictional drag, Pig. 6*5 whereas the coefficient of friction 
decreases significantly, Fig. 6 •4* Ihe percentage rate of increase 
or decrease is higher for higher values of N. In case of longitudinal 
routines s there is a slight decrease in load carrying capacity, 

Pig. 6.2, and a slight increase in fidctional drag, Pig. 6.5 
consequently the coefficient of friction increases for all values 
of N. Again it is seen that the rate of increase is further enhanced 
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Fiy. 6-5 Load parameter vs. length ratio j.arameter L for 
various values of coupling numter N(B=0'49, 
longitudinal roughness) 




LONGITUDINAL ROUGHNESS 



Frictional drag parameter Ff, vs. length ri 
L tor various values of coupling number N 
lonuitudinal roughness) 



LONGITUDINAL ROUGHNESS 



Coefficient of friction parameter vs. length rat 
parameter L for various values of coupling numiber 
( B =0*4 9; longitudinal roughness) 
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Load parameter \A/r vs. length ratio parameter L for various values of 
number N |B = 0*49, transverse roughness) 


TRANSVERSE ROUGHNES 
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F'ig.6-9 Frictional drag parameter Ff:, vs length ratio parameter L 

for various values of coupling number N (B:0-49 transverse 
roughness ) 
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B33‘ ^ 24 36 48 L - ^ 60 72 84 96 

pi^p.-lO Coefficient of friction parameter Cr vs. length ratio parameter L for various value 
of coupling number N (B = 0'49 transverse roughness) 
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CHiPTEEl ¥II 

SQUEEZE Ein«B BEW7ESH ROUGH SURFACES 

7.1 i:iDROUUC?IO!T 

A nuEiber of experiments in which the evidence of rheological 
ahnormalities, like enhancement of viscosity, has "been rather 
extensive and convincing, are of squeeze film type. For example 
Heeds [ 22] conducted a series of experiments to detect the 
influence of boundary surface on the viscosity of lubricant films 
between two optically plane parallel circ^^lar disks approaching 
each other with no tangential velocity. He measurod film thickness 
down to 0.000655 mm. For values less than 0,00127 mm an increasing 
discrepancy appeared between measured and theoretical (as predicted 
by the classical Hewtonian theory) intervals of time. It actually 
took longer time for plates to approach than what was predicted 
by the theoiry, indicating that the effective viscosity of the 
lubricant in such thin films had somehow increased. Effective 
viscosity in the boundary film nearly five times the bulk value was 
observed. Ho conclusive explanation could be offered except to 
suggest that the proximity of the metal surfaces influenced the 
viscosity of the fluid, causing the fluid to become more rigid. 

Experiments by Puks [28,29] also fall under the same 
category. Various fluids (like solution of pure fatty acids in 
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ptire hydrocarbon solvent) were sq.ueezed between two plates 
submerged in the fluid and measuring the gap bet'.'aen the plates as a 
function of time. It was observed that even after squeezing for 
a very long time a residual film remained between the plates. 

I'Jora recent experiments of this type are by Drauglis ot.alL 31 ] 
and Askwith, Cameron and Gohar [ 36] . In the experiment of Drauglis 
et,al, results similar to Fuks [ 28,29 ] were obtained, residual 
films were observed only in the presence of additives. However in 

Askwith et ,al [ 36 ] squeeze film studies showed that the plates 

-4 

came to rest at a separation of about 2 x 10 mm with all fluids, 
even pure Cetane. 

It has been suggested that aspeidties on the solid surface 
might account for the residual film . Also dirt and undissolved 
materials in the films have been invoked as the explanation for 
this residual film [ 30 ] • According to Drauglis et .al [ 51 ] 

"if the data of Fuks are all to be explained in this ms^nner, then 
one must postulate dirt which changes its properties vrlth. temperature 
length of molecule of additive and concentration of additive but such 
dirt would be a most -unusual material" . 

A proper theoretical explanation of the residual film reported 
by various workers is therefore clearly needed for a better understanding 
of thin film lubrication technology , especially vfhen the height of the 
■ roughness asperities is of the same order as 'that of the thin films . 
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Ifotivated hy these considerations the squeeze film 
problems have been studied by various authors [ 55j 68-71j75»76 J 
from aicroctmtinuum view point. Inherent in all earlier theoretical 
analyses for squeeze films including those from the microcontinuum 
view point v/as a fundamental assumption of perfectly smooth surfaces. 

It has long been recognized that owing to machining limitations, 
this assumption is rather unrealistic, particularly in bearings 
where the bearing clearance is very small. Of interest, in this 
direction are the works of Prakash and londer [ IO4] who studied 
squeeze films between rough circular plates and Prakash and 
Christensen [ IO5 ] who studied squeeze films between rough 
rectangular plates. In analysing the rou^ness effects use was 
made of the stochastic models for hydrodynamic lubrication of roxi^ 
surfaces developed by Christensen and Tender [ 99>100] and Christensen 
[ 101] for Fewtonian fluids. 

Hence in this Chapter, the lubrication theory for mioropolar 
fluid for rough surfaces, advanced in Chapter 71, fs applied 

to the squeezing between parallel plates. Three standard geometries 
have been analyzed, i,e, two dimensional parallel plates (Pig,7,l(a))j 
circular plates and rectangular plates (Pig,7«l(b)) (three dimensional) . 
It is our strong belief that this theory could provide a proper 
theoretical explanation to the existence of residual films. 

The application of this theory may be found in the lubrication 
of rough surfaces where long chain length additives are used in 
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lubricants. An example to this context is the slrndding of a rolling 
pneumatic tyre, even at low vehicle velocities on a wot surfsx;e. 

The main feature contributing to this slddding is the existence of 
a squeeze film in the contact area between tyre and surface. This 
could be done by minimizing the sinkage time. The beha-viaur of a 
tyre element in this zone my be idealized by the psra,llel plate 
squeeze film onto a rough airface and the fluid in between as a 
micropolar fluid on account of suspension of rigid particle additives 
(e.g. dirt particles) . A theoretical and experimental study of such 
a situation has been made by Moore [ 106 ] , where the asperity height 
of a randomly rough surface have all been assumed same. 

Further application of this study may be found in human 
joint lubrication where the mechanism is of squeeze film type. 

The lubricating synovial fluid in joints can be considered as a 
micropolar fluid and the cartilage airfaces rough ( a detailed 
justification is given in Chapter VIIl) . 


7.2 MiTIISI'aTICAL MALYSIS 

The generalized forms of Reynolds equation for rough surfaces 
from the micro continuum approach have been derived in Chapter VI, 

For squeezing, equations (6 .11), (6,15) s.nd (6.19) reduce, respectively, 
to the following , 


Longitudinal, one-dimensional roughness, 
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3 

ax 


at 


Transverse, one -dimensional rougiiness, 


{(E(f(H,£,hr^))“^-~^} {E(f(lT,£,h)) 


as 






BE(h) 

at 


Isotropic roughness, 


{E(f(N,£,h)) {E(f(N,£,h)) 


*p- 


az 

aE(h) 

at 


where 


h = h^(x,z,t) + \(x,z). 


(7.1) 


(7.2) 


(7.5) 


E( )= / ( ) f(h^) dhg 


f(h ) is the rougiiness distrihution function of the random 
s 

variable h given by eqn. (6 .20), Bf and £ are defined in (2,18) 
s 

and f(N,£,h.) by eq.n. (2 .55) • 

Since the roughness distribution function is assuned 
to be symmetric, 


(7.4) 

(7.5) 


E(h) = h. 


n 


(7.6) 
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7.5 SqUEBZE FiaS BEPWEEII ROUGH I1TPII7CTELY LOITG PAEALLBL FIAIES 

(i) Longitudinal < One-Dimensional Rou^ness 

Eor this type of roughness, the ridges and valleys are asstmed 
to he running along x-direotion. 

The film thickness assumes the form 

D- = \(t) + hg(z) ( 7 . 7 ) 

and the one-dimensional Reynolds equation, from eqn, (7 »l) , is 

^ ( 7 .B) 

The pressure boundary conditions are 

E(p) = o, at X = ( 7 * 9 ) 

= o (e(p) maximum) , at x = o (7 .10) 


where Is the half x-dimension of the plate. The mean pressure 
distribution is 


dh 


_/ \ 1 .. ““n / 2 ^ 2s 1 

E(p) = - 2 IT E(f(R,il,h5T 


(7.11) 


The expected value of load carrying capacity per unit width is 

£1 

E(W) = 2 / E(p) dx (7.12) 

o 


or 


dh 


E(W) = 


51 


n 


dt E(f^(i,£,h)) 


(7.15) 
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The non-dimensional load carrying capacity is 


E (?0 = 


bC>7) 


P; dh 

/ 1 Un 

^■“3 v5 dt ^ 




( 17 , 1 , 3 ) 


( 7 . 14 ) 


where 


3 


P^(1T,L,B) 


323' _B 


(1+HJ5 (1+h ) 


12 


ooth S(l+E^) } 


( 7 . 15 ) 


and 


h h , 

H = ^ , L = -ii , 3 = f 

s \ \ 


(7.16) 


The time of approach, t, from an initial (nominal) film thickness 
h^ to a final (nominal) film thickness 1^, is given hy 




0 3 h 

„ y 2-:; n 
2 1 j 


dh 


n 


3 e(W)J ”“B(f (iTjJl ,h) ) 

^i 

The non-dimensional time of approach is 

H SB. 


T = 


1 ( 2 _ l.A 
^2^ 3 E(W) ■' 


r n n 

i 


( 7 . 17 ) 


(7.18) 


where 
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B _ (h^+hj5 (h^+hJ 

■iV / < 


32 B' _B 


n s'^ + ' 

*12 




ooth (H « ) ) dE 

9 ^ Y\ a' c 


2L 


(7.19) 


and 


h h ^n. 

Tj ^ n -q- _ s Y a ^ B 

n ’" h ’ ■"'s h ’ ~ i * ~ h 

ri. n. n. 

IX X 


(7 . 20 ) 


(ii) Transverse, One-Bimensional Roughness 

In this roughness, the ridges and valleys are assumed to 
be running perpendicular to x-direction. 

The corresponding Reynolds e<luation, from eq.n»(7.2), 
is 




dh 


(7.21) 


where 

h = hj^(t) + hg(x) 

Proceeding in a way similar to longitudinal rouglmess 
the average value of squeeze film characteristics are given 
by 


o t: 

E(w) 


-) 






- 3 WWJ ^^f(IT,£,h) ^ "^n 


( 7 . 22 ) 


(7.23) 


(7.24) 
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The dimensionless load capacity and squeeze time are 


eM = 


H 

T = - r dH^ 


(7.25) 

( 7 . 26 ) 


where 


^ dH 

P,(N,L,b) = ^ ^ 

^ 8B^ _3 , 12(1+H ) 

{ (l+H )^+ 


6N(l+H y _ 

jj-^oth~(l+Hp > 


(7.27) 


(B^-H^)^ dH 


P4(H,L,H^,B) = 


B 

105 j - 

, 6H(H +H )^" 

{ (VH^)hI2(H^^s^)--..■ I 3, rf:hS(H^4B^)} 


( 7 . 28 ) 


and the other non-dimensional parameters for load capacity are 
defined according to eqn. (7.16) and for squeeze time as per 
eqn . (7 .20) , 

(iii) Isotronio Roughness 


The roughness is uniformly distributed over the surface 


of the plate. Bor two dimensional geometries the isotropic 
roughness will he equivalent to longitudinal roughness. 
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7.4 SQUEEZE Fim Bm'mmi rouCtH circtjlar piates 

(i) Radial, One-J>iiaensiorLal Roughness 

lii this model, the roughness is assumed to have the form 
of long narrow ridges and valleys running in radial direction, and, 
therefore, film thickness assumes the form 

h = h^(t) + h^(e) 

!Phe governing differential equation, from eqn, (7*1) in 
polar co-ordinates, is 


i|^{r E(f(U,£,h)) 


1 3 
“2 86 


9h 


(E (f(N,Ji,hrt 


For an axisymmetry case, eqn. (7 »30) reduces to 


dh 


•i^CrE (f{H,i,h)) 


The boundary conditions for average pressure are 
E(p) = o, at r = r^ 

= o (e(p) maximum ) , at r = o 


where r^ is the radius of the circular plate. 


(7 .29) 


(7.50) 


(7.51) 


(7 .32) 


(7.35) 
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The mean pressure distrihution is given hy 


E(p) 


IL 

dt ^1 ‘ 


t 


1 

E(f(]j,£,h)y 


(7.34) 


The average instantaneous load capacity and time of syi^ezing 
are 


](¥) = J 2ir r E(p) dr 
o 

ydh 


JL " 4 1 

8 dt 1 ~'E(f(N,£ ,h) ) 


t 


4 ^ 

V y r^ 3 

- iE(fr h 


dh 

n 

E(f(N/,h)) 


The non-dimensional load carrying capacity and squeeze time 
are 


(7.35) 


(7.36) 


bTwT = 


E(¥) 


dh r!; 

/ 2L!i E -i'l 

8 dt .y 

h 

n 


1 

F^(1T,L,B) 


t ^n 

^ 8 KfJ ^”2 


(7.37) 


(7.58) 


where P^(n,L,B) is defined hy eqn, (7*15) and E2(N>L»H^»B) hy 
eq.n, (7.19) and the non-dimensional paraneters for load are 
according to (7 .16) and for squeeze time as per eq.n. (7 ,20) , 
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(ii) Circumferential, One-Simensioml Roughness 

In this model, the ridges and valleys are assumed to he 
running along 6 direction. The film thickness is 

h = h (t) + h (r) (7.39) 

iX o 

The corresponding Reynolds eq.uation, from eqn. (7.2), 
in polar co-ordinates, is 

+ -2-5r <E{f(v,h)) -fish} (7.40) 

X 

For an axisymmetry case eq.n. (7 .40) reduces to. 


{r(E((f(R,£,hr^ ) 


= y 



(7.41) 


Proceeding in a similar way as in the case of radial roughness, 
the load capacity, squeeze time and their non-dimensional forms 
are given hy 

_ ydh 

E(7r) = - g "JF ^1 ®^f (ii,£ ,hy^ 

4 r 1 

^ “ ai(f) ^1 ( 

n 

X 


) 


(7 .42) 


(7.43) 
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and 


E(W) 




( - 


n 


P^(1T,L,B) 


H 


T = 


n 


Trp r:: 


= - / F,(H,L,H ,B) dH 
J 4'' » » n 


8h S(W) 
i 


where Pj(N,L,B) is defined "by eqn. (? •2?) and (]:T,L,H^,B) 
hy eqn. (7 ♦28) and the other non-dimensional parameters for 
loa^ capacity are defined hy eqn. (7 •16) and for time height 
relation hy eqn, (7#20) . 


(iii) Isotropic Roughness 

The modified form of Reynolds equation, from eqn, (7 #5) » in 
polar form, is 

i|^{rE(f(H,«,h)) tE(r(:r/,h)) 


3r 


236 


38 


3h 


= V 


n 


3t 


For axisymmetry case the form of isotropic rou^ness will 
he identical to radial roughness. 

It is evident that the qualitative and quantitative 
effects of roughness will he same for both the geometries, i.e. 


(7.44) 


(7.45) 


(7.46) 
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inrinitely long; parallel plates and circular plates (axisyimnetry 
case) • Hence the results are given only for circular plates. 

7.5 SQ,UEB2;iITG BmPMEEB ROUGH EEgTAmTJLiR PIATES 

The previous sections were concerned with the one-din ensional 
prohlems. Therefore the side leakage effects, the width/length 
ratio effects for unidirectional roughness were not detected. 

Moreover the distinction between the isotropic raighness and 
the one-dimensional roughness could not be made. Hence, in this 
section, the problem of sq.ueezing between rectangular plates is 
analyzed in order to stui^ the effects of aspect ratio (width/length 
ratio) on the lubrication of rough surfaces with micropolar fluids. 

The geometry of the system is shown in Fig.7»l(b), 

Since there is no relative tangential velocity between 
the surfaces and because of symmetry, the one-form of one-dimensional 
roughness can be obtained from the other, simply by a rotation of 
co-ordinate axes. Hence, only the case of one -dimensional 
roughness, with roughness having forms of long narrow ridges 
and valleys running in the x-direction is studied. This will 
be termed as 'Unidirectional Roughness’ , The case of isotropic rou^ness 
is also analjTaed, 

(i) Unidirectional Roughness 

Since, in this model, the rou^ness is assumed to have the 
form of long, narrow ridges and valleys running in the x-direction. 
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the film thiclcness assumes the form 

h = h^(t) + hg(z) (7.47) 

The corresponding Reynolds type equation, from eq.n. (7«l)j is 


ax 




3h 


= y 


n 

a t 


(7.48) 


Equation (7.48) can be written as 


1^S(p1 , .2 _3.^.E(^ . 


ah- 


n 


ax 


R 


a z 


E(f(R^ ,h))' at 


(7.49) 


where 


“^R “ “R 


6^(11, Jl,hj = {E(f(N,Jl,h)) X y ) 


1 

‘2 


(7.50) 


Equation (7*49) is Poisson-type equation and can be solved usual 
separation of variable techniques with the boundary conditions 


E(p) 

E(p) 


o at X = 0 and x = 

O at 25 » + 


(7.51) 

(7.52) 


to give the mean pressure distribution 


E(p) = 


. 2 dh 

4 V ap n 

dt 


cosh 


— I 


nTT z 
^2^R 


cosh 


nw b. 


- l)x 


srn 


nir X 


(7.55) 
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where 8-2 is the x-diaens ion. and h2 4s the z-dimension of tlie 
rectangnlar plates. 


The man instantaneous load carrying capacity is 
a 

2 

E(¥f) - j j S(p) dx dz 


or 


I6y dh 

E{W) - -J- (-^) 


(. — a tanh 


nTT K 


1. 


n=l,5,5 n^ir 


The non-dimensional load carrying capacity is 
E(¥) IT 4^5 


W) = 


16 VI A' 


n 


o dh 
2 n 


dt 


(7.54) 


(7.55) 



1 

P^(N,L,B) 


“ 6,3 nirv, V 

n=l, 5}5 


(7.56) 


where 


^ E “ f^:P^ » 4. - ^ 

P^(H,L,B) is defined hy egn. (7.15) » ^’^(l^jL,®) ly eqn. (7.27) 
and other non-dimensional parameters, i*e. B,L, are defined 


(7 .57) 
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according to eq.n. (7»16) . 
The time of approach is 


48 U a-r 


h 


n 




/ 


h^ S(i(TJ,£,h)) 


I 


/ 

( — — ~ tanh 


ntr "b^ 


n=l,3>5 ’n'li"' 


The non-dimensional time of approach is 


T = 


t TT £(?) h^ 


48 y A 




2h 


4^ 


(7.58) 


+ v"^ / 


A H 
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nirv 

— ^tanh — 12: 

n=l,3,5 a 26 


I ( 


A 


R, 


2n 


.)dH 

4 n 


■where 


—2 

R 



(7.59) 


( 7 . 60 ) 


The paramebers H ,B,H and L are defined according to eq.n. (7 .2o) 
n s 

and F 2 ( 1 ^jL,\,B) by eqn. ( 7 . 19 ) and P^(R,L,H^,B) by ©in. (7.28). 
(ii) Isotropic Roughness 


The film thickness is 


h = \(t) + hg(x,z,e) 


(7.61) 
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and the Reynolds-type etliiation for this type of rotighness is 




}+ — {B(f(lJ,£,h)) 


9 z 


9\ 

~W 


(7.62) 


or 


3. E(p) 4 ^. 5(p) ^ 


9X 


9z 


y n 

E(f(H,j,h)) dt 


(T.S« 


The various bearing characteristics and their non-dinensional 
forms are given as 

Pressure distribution, 


4y a^ dh. 


cosh 


n'lT z 

3#r 


" ^5 dt"" n-l,3,5n3^ 


tlTT X 


oosh 


2a^ 


(7.64) 


Load capacity, 


dh 


^ dt E(f(N,"£,h)) 


(7.65) 


where 


16 y a^ “ a nvb^ 

K = -7-^ I ( -~ tanh - 2^ 
TT^ n=l,5,5 rP 2 


2n^ 


( 7 . 66 ) 
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and is independent of the rou^ness parameter. 


Squeeze time, 


t = 


K 


n dh 

I 


h S(f(l,£,h)) 

n. 

1 


(7.67) 


The load capacity and squeezing time can be put in non-dimensional 
fomn as 


S(¥) = 


E(W) 


dh , 
dt n 


FiC¥,l,b) 


(7.68) 


T = 


K 


— / 


n 


dH 


n 


E(W)h' 




(7.69) 


n- 

X 


Thus it is evident that the qualitative effects of isotropic 
roughness on rectangular plate squeeze films are identical to that 
of radial one - dimensional roughness in circular plates (or 
longitudinal roughness in infinitely long parallel plates) * 


7,6, liDIviBRlGAL RESULTS 

The roughness effects on various bearing characteristics 
are studied through the following parameters 
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(i) Load Capacity Ratio Parameter 

The load ratio parameter, is defined t; 


e07)- (s(77)).. 


=0 


E 


dlwj) 


B=o 


(7 -70) 


where CE(W))g_^ is the load capacity for smooth sixrfaces with the 
same nominal film thickness. 


(ii) Time Ratio Parameter 

The time ratio parameter, T^^, is defined by 


■R 




(7.71) 


Where (T)b_q 3:epresents the squeeze time for smooth surfaces 
to reach to the same nominal film thickness. While calculating 
this ratio the load is assumed to he same for rough sarfaces 
and smooth sijirfaces as well. 


7.7 RESUITS ML LlSCUSSIOh 
(i) dimensionless Parametacs 

In addition to the well known pairameters for smooth squeeze 
films , lubricated with Newtonian fluids, four new paranffiters 
N, 1, B and are introduced in the present study. 

The mioropolar parameters N and L have been discussed in 
Chapter III and the roughness parameter B in Chapter 71, Here it 
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suffices to say that the effects of surface roughness Trould he 
pronounced for large values of B and N and snail values of L, 

^2 

The non-dimensional parameter is termed as aspect 

Jx 0^2 

ratio. In case of unidirectional roughness the rolo of this 
parameter might be impoirbant because the value of will determine 
uhether the nature of the unidirectional roughness is longitudinal or 
transverse. 

(ii) Bearing Characteristics 

Bigs, 7*2 to 7*4 3.xe the graphs for circular plates 
sq.ueeze film and Pigs, 7*5 to 7*7 rectangular plates squeeze 
film. 

Further, it should be noted that the cases of radial 
and circumferential roughness for circular plates are equivalent 
to the cases of longitudinal and transverse roughness, respectively, 
for the infinitely long parallel plates. 

Fig, 7*2 is the graph of vs, B where the nominal 
film height is kept constant. It is seen that when B o,5 
the load capacity for rough surfaces does not differ significantly 
as compared to smooth ones. Even for micropolar fluid this 
assertion is tmo. Beyond B >■ o*5 the rou^ness effects become 
prominent and as B ^ 1 the load capacity is affected significantly. 
It is seen that the radial roughness decreases the load capacity 
while 0 incumfer ent aal rou^ness increases it, The load capacity 
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for circumferential roughness may he times and. for radial roughness 

O 

times the load for corresponding smooth surfac<=s, as B 1, 

These limits are for Ne'"tonian fluid and can he obtained directly 
from eqn.s (7.44) and (7.37) t respectively hy taking the limits 
B -i- 1 and IT o. For micropolar fluids the rate of increase or 
decrease in load capacity is higher and the deviations from Newtonian 
rate are more prominoit at higher values of rou^ness parameter B and 
higher values of N» 

Fig, 7,3 and Fig. 7,4 are the graphs for Tg^ vs. non-dimensional 
film height for B = o,l and B = o,4j respectively. From Fig. 7*5 
it is clear that the effects of roughness are prominent only when 
the roughness paramsbec is of the same order as the film thickness, 

T} 

i.e, when ^ approaches unity. It is also evident that the circumfer - 
n 

ential roughness increases the sihkage time while the radial roughness 
decreases it aid the effects of circumferential roughness are more 
prominent than the radial roughness* The reason hoing same as that 
given for load capacity. 

Fig, 7,4 shows a similar trend as shown hy Fig, 7,3, except 
for the accentuation of rou^ness effects. 

Fig, 7,5 is the graph between and B, It is seen that 
the roughness effects are significant only for B > 0,4 and for 
values of not around 1,5* 5’° 3? values of around 1,5 the 

roughness effects are sognifioant only for B very close to unity. 
Further, it is seen that the load capacity for rough plates is higher 
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than that for smooth plates for smaller values of Vj^ (i.e. v < I. 5 ) . 
For values of from v = 1,5 to ~ 2.o, approximate!”, the load 
capacity decrea.ses first, reaches a miniraum, and then increases. 

For values of higher than 2 ,0 the load capacity decreases for 
almost all values of 1, A mathematical reason for ITewtonian fluids 
is given in the following (for micropolar fluids a similar reasoning 
can be advanced, however, the mthanatics involved will be complicated 
and hence only numerical results are given in this case to illustrate 
a similar phenomenon) . 

She unidirectional rou^ness pattern, with roughness ridges 
and valleys running in x-direction, tends to increase the x-component 
of the flow, E(q_^) , whereas the flow in z-direction, E(q^) , is 
reduced, mathematically, 






3 n 


(7.72) 


and 


B(q^) 


1 32b^ 

lai 35 


[ 3(5bW)(b^-h^) log 


h +b 
n 




n 


+ 2 bh^(l 5 hJ- 13 b^)]“^ ( 7 . 75 ) 

It is obvious from eq.n . (7 *72) that the flux in x-direction 
increases and reaches to a maxiraum value as b -»■ h^^ (or B ->■ l) , which 
is times the corresponding flux for smooth surfaces . The flux 
in z-direotion decreases as~^( « B) increases and as b h^(or B -v l) 
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3 

it becomes — ^ times the corresponding flux for smooth surfaces. 

Further, it must he noted that for a constant area hearing 
the component of the flow along the smaller dimension would 
he larger because most of the fluid would have a tendency to escape 
along smaller dimaision. Thus if ^ Sc ^ smooth surfaces) , 

Thus when v a> v (v ^ 1*5» v sH-y he termed as critical 
aspect ratio ) hut not more than 2 ,o (approximately) , the x-dimension 
is smaller than the z-dimension and the increase in more 

than enough to offset the decrease in E(q^) , This remains true 
for smaller values of B, For higher ■values of B, on the o’Kher 
hand, the increase in will not he able to balance the 

decrease in E(q. ) . However when v. is sufficiently large, 
i»e, more than 2,o (approximately), the effect of decrease in 
E(q^) will not he more than the effect of increase in E(q^) , for 
a constant area hearing, 

Also, the rou^ness effects are more pronou.nced for 
micropolar fluids, i.o. the rate of increase (for smaller values 
of v^) and the rate of decrease (for higher values of v^) are 
hi^er when mioropolarity of the fluid is taken into consideration. 

Fig, 7,6 is the graph, T^ vs. nominal film height 
keeping the roughness parameter fixed (B = o.2). It is seen that 
for sfliialler values of the uni-directional roughness increases ihe 
sinkage time, i.e . the behaviour is similar to that of transverse 
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roughness. It is seen that around = 1«5 the roughness effects 
are not prominent . ?or considerably high values of v^, the 
uni-directional roughness decreases the sihkage tine for all values 
of B . 

Fig. 7*7 is the time - height relation for B = 0,4. 

The trend is similar to Fig, 7*6 except for the accentuation of 
the roughness effects due to the micropolarity. 

S'or isotropic roughness the load capacity and sinkage 
time are always lower than the c ctrresponding smooth ones. It is 
seen that after ^ I 5 , the unidirectional rou^ness and isotropic 
roughness yield approximately the same results. The case of isotropic 
roughness from uni-directional rou^iness can be derived by letting 

It should, however, be remembered that the comparisons 
here have been made on the basis of nominal film hei^t . Therefore, 
if h^ is the minimum film thickness for smooth surfaces, 
would be minimum film thickness for rough surfaces and hence the 
minimum hei^t permissible in order to avoid a surface - to - surface 
contact would be much less than that in the case of smooth ones. 
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CHAPTER VIII 

LUBRICATION OP HIRviAlJ JOINTS 


8.1 IIITRODUCTICN 

Two centuries of study has failed to unravel completely 
the mysteries of the mechanism of human joints, although from the 
engineers point of view they represent one of the closest points 
of approach "between the inventions of man and the work of nature. 

The synovial joints provided "by the nature in the human 
body to carry out the trouble free motion of one bone past another, 
have long been identified as bearing systems. With the exception 
of the joints of the middle ear and that between the lower jaw and 
the skull, all synovial joints are potentially weight-bearing, 
although some exercise this function more often than others. These 
joints function as excellent bearings in the biological conditions 
and generally support considerable load, while providing, low friction 
service over a long span of one^ s life time. There, exists no mechanical 
device which yields as low a coefficient of friction as that found 
in the sjnaovial joints. Thus a clear visualization of the synovial 
joints is likely not only to provide a guide-lino to improve the 
existing mechanical "bearings but also help in the replacement of a 
faulty hximan joint. 


The human joint, of course, is an excessively complicated 
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system from an engineering viev.” point and it is not possiMo to 
build as a mechanical replica in which all the features of the r»l 
thing are represented, nor indeed is it possible to dOTice a 
mathematical description that is conpletoly comprohonsivo . It is 
more practical and probably more logical to oommonco inTestigation 
in the opposite manner, i.e. to start with the most simple squeeze 
film and to introduce complicating features one at a time, until 
either a situation of sufficioat complexity is reached or the engineering 
and computational aspects inhibit further progress* 

This fundamental approach has led to the development of 
variotis theories of human joint lubrication. Serious study of the 
lubrication mechanism in human joints dates back to almost half a 
century ago. MacGonaill [ IO 7 ] suggested as early as 1932, that 
the synovial joints function as the hydrodynamically lubricated 
bearings, Chamey [ 108 ] challenged this view in 1959* Sis 
experimaatal evidence supported a boundary lubrication action, 

Chamey observed a chemical affinity of hyaluronic acid molecules 
with the cartilage surface, McCutohen [ 109-112 1 introduced the 
porous nature of the cartilage into the lubrication concept, when 
he pm posed his "weeping lubrication" ideas. Dint enf ass [ 115 ] 
put forward the concept of elastohydrodynamic lubrication for 
synovial joints, 'sdiich was later supported by Tanner [ II 4 ] « Any 
possibility of boundary lubrication was ruled out by Dintenfass, 
while Tanner indicated a sli^t possibility of boundary lubrication 
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at the asperity contacts . 

However, it was only in I967 that Dowson [ 115 ] realized 
that one mode of luhrication is insafficient to explain the joint 
laeclmnism. He concluded that "the major l\ibrication mechanism 
TiTould seem to he some fom of elastohydrod^/namio action determined 
by sliding or squeeze filnwiction between porous surfaces with 
boundary lubrication providing the surface protection in cases of 
seTere loading and little movemait" , 1 little later in I968, 

Walker et.al [ 116 ], suggested a new concept of joint lubrication,!* e . 
"Boosted lubrication". According to this concept of lubrication, 
squeeze-film action leads to a concentration of hyaluronic acid- 
protein complex in the lubricant as a result of diffusion of water 
and low molecular weight substances through the porous cartilage, 
and the restricated gap between the approaching cartilage surfaces. 

The increased concentration of hyaluronic acid molecule would 
give rise to an increase in viscosity of the synovial fluid in 
accordance with the findings of Negami [ II7 ] and the proposal 
is also consistent with the formation of gels on the cahtilago 
surface (Maroudas [ 118 ] ) . Other work in this direction is by 
Higginson and Norman [ 119] ''i'^h.o inferred that elastohydrodynamic 

luhcication plays a significant role. IJnswo3d;h ot ,al [ 120] and 
Wright and Dowson [ 121 ] have strongly suggested boundary lubrication. 

Though, Dintenfass [ 11 , Tanner [ II4 ] and Dowson [ II5] 
had made some numerical calculations for the fluid film thickness, 
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Fein [ 122 ] was proba-bly the first to give a mathc-Entical Eodel for 
synovial joints • The joints system was considered to he elastically 
eq.nivalent to sphere-on-a flat system. Mow and Ling [ 125 ] made 
use of Biot’s formulation of poroelastic theory for giving the 
displacement field equation and also the strain displacement 
relationship. The fluid flow in the porous matrix vjas taken to he 
governed by simple Darcy's law, Dowson et ,al L 12/1-125 ] gave yet 
another simple formulation and accounted for viscosity variation 
with the hyaluronate concentration. The boosted squeeze time was 
found to bo in good agreement ^7ith the experimental result. Ling 
[ 126 ] presented an. analytical model to account for the non-linearity 
of the cartilage. It was pointed out that the two mechanisms of 
'weeping', and 'boosted' are not exclusive. Higginson and lorman [l27 ] , 
Mansour and Mow [ 128,129 ] and Higginson [ 150 ] have given more 
rigorous mathematical models. They considered the elastic as well 
as the porous nature of the cartilage. 

It is however observed that the present situation regarding 
joint lubrication is not very clear, Sarlier works based on the 
classical continuum theory, could not explain the observed effective 
increase in viscosity near the solid boundary. 

Little is known about the rheological nature and quantity 
of lubricant present in the human joints. Such knov/lsdge is 
■Q^^i^gQUtly required for the development of materials and desxgns and 
the better understanding of the tribological characteristics of 
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replacement hearings for thehuran body. 

It Vk’ould not be out of place to mention a little about 
the synovial fluid. Normal synovial fluid is generally clear or 
yellowish and viscous. It can be briefly described as a dialysate 
of blood plasma. It contains about one third of the protein 
concentration of the plasma. In addition the synovial fluid contains 
a very important polymer known as hyaluronic acid or sometimes called 
mucopolysaccharide. It is this hyaluronic acid molecule #iich gives 
synovial fluid its slippery and stringy behaviour. It also gives 
the characteristic non-Newtonian behaviour [ 15 I ] • It was concluded 

by McCutchen [ 111 ] that the synovial fluid is a boundary lubricant , 

The synovial mucin are polymers of high molecular weight. Estimates 
of the molecular weight of the ^movial mucin are as high as several 
million. Chemically, it is hyaluronic acid conibined more or less 
firmly with less or more protein. 

Thus in studying the human joint lubrication problan, one 
must examine not only the continuum behaviour of the synovial 
fluid but also the microscopic structure. In obtaining the constitutive 
equations one must, therefore, resort to the miorooontinuum approach 
as the continuum approach will be less definite about the detail 
structure of liie material , 

This provides the motivation for the present work. The 
micropolar fluid model proposed by Eringen [44} 


t 


serves as a 
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satisfactory model for the description of the flovr hehaviour of 
the polymer fluids and roal fluid suspensions. It h^s alroady 
heen pointed out that the micrcpolar fluid theoiy can be utilized 
to explain the rheological behaviour (increase in viscosity and 
boundary lubrication) of the lubricant. It also servos as a 
satisfactory model, uhen the fluid has long chain additivo 
molecules , 

It is, therefore, noted that the synovial fluid may be 
thou^t of as a micropolar fluid because of the long drain hyaluronic acid 
molecules. The increase in the concentration of the hyaluronate 
in the synovial fluid, which is the important feature of the boosted 
lubrication, causing an increase in the effective viscosity, can 
also be accounted for in this theory, by way of increase in the 
mioromotions and the couple stress. In this analysis, it is 
assumed that the human joint can be approximated by a system of a 
porous spherical bearing, Micropolar fluid theory is applied. 

Slip which almost invariably occurs near the porous boundary is 
considered by using the Beavers and Joseph [ 152] velocity slip 
criterion. The elasticity of the cart ilage ,i .e , the porous matrix, 
is ignored to avoid the additional complexity in the already sufficiently 
c complex model. Thus the present analysis may be termed as a first 
approximation to the hiomn joint lubrication theory, using the 
miorocontinuum approach. The only work in this direction is by 
Tondon and Jaggi [ 75»74 1 • have used, rather too simple a 
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geometry for the human joint. The validity of various equations is 
dehatahle as vi-ill he obvious from the present analysis. The 
coeffic ient of friction is calculated vrongly. Moreover the full 
potsatial of the micropolar fluid theory has not been realised, 

8.2 MATOMT OF A SIMOYIAL JOINT 

The different foxjns of synovial joints in the Human body 
can not be readily reduced to a single model. For e3cample, tbe 
hip joint is a close fitting ball and socket affair, to be 
appcoximated by an equivalent spherical bearing, while the ankle 
joint can be approximated by a cylinder-cylinder arrangement, 
resembling a close fitting journal bearing. Fig, 8,1, shows a 
joint illustrating the saliait features encountered in synovial 
joints. In this schematic repre saltation of a human joint, the 
load transmitting structural menbers are the bones. The ends of 
the bones ‘vdiich are usually globular in appearance, are covered with 
a layer of relatively soft and porous material called cartilage. 

The articular cavity contains the sjmovial fluid, which is 
responsible for the coii5)lex joint lubrication phenomenon, 

8.5 THEOBETICAL FOEMDIATICfr „ ASSDICPTI0H5 AID EQUATIONS 

Geometrically the system of synovial joint, especially 
the hip joint may be represented by a sphere on a porous surface. 
Both the sphere and porous surface are assumed to he rigid. The 
simplicity of the arrangement is shown in Fig, 8.2, The noimal 
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Fig 8.1(a) Synovial joint 
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Fig 8.Ub) Hip joint 
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approach of the sphere is achierved hy allowing it to fall under 
gravity on to the lower porous surface. The syc’:er. is assumed 
to he lubricated by a highly non-lTewtonian *Syncvial Fluid* which 
has long chain hyaluronate molecules . Micro polar fluid theory is 
assumed to govern the behaviour of the synovial fluid. 


(i) Sasic Equations 


The basic governing eqriations for a three dimensional 
flow field under the assumptions stated in Chapter II, from eqns, 
(2.52) to (2.55) ,are 


(y+ fx) 

9 y 


3 v„ 
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3 x 


= o 


( 8 . 1 ) 


(u+ Jx) ^2“^' 
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= o 


( 8 . 2 ) 


a^Vq g 

^2Xv 

3/ " 


= o 


(8.3) 
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The boundary conditions are 


u = TJ^^, ^ ^ = '^5 “ o, at y = o 

/T 


(8.5) 


u = IJ, 
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y«h » ^ “ ^22 
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y=ii 


= o, = o, at y “ h 


( 8 . 6 ) 
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where (|)„ is the porosity of the cartilage surface and A is the 

Cu 

slip paraneter. Condition (8.5) is the no-slip boundary condition 

for the velocity and Fdcrorotation velocity at the non-porous 

boundary and condition (8.6) is the slip boundary condition for 

to 

velocity according Beavers and Joseph [ 132 ] , however, the 

microrotation is assumed to vani^ at the porous surface also . 

(ii) Velocity distribution 

Solving equations (s .l) to (8.4) velocity components 
u,w and microrctation velocity components ; with the boundary 

conditions (8.5) and (8.6)> we obtain 

^-^2 I? “ \iy^ - ? ®iiy ^ ?^^®3i ®4i ®2i 

(8 .7) 

_ 2 0 p - 2 ^2 

» - irC-f ys - ~ - V (*32 stah “y + V ““sii "y)+h2 

(8 . 8 ) 

^1 “ tI ” \ 2 } “ ®52 ” ®42 

^3 “ " 2T^y ‘ax " ®ii^ ®3l ®4l ^ (8.10) 

where 

V [ (n^,-n. 2 ) + I: 2* ^ + h - f-^tarih f) ] 

{4^(l-ll^) + h - ■^tanh 


(8 .11) 
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-5 p 

r-rr — r — {B (cosh sih-l) + h -r — / (8.12) 

2i xl pm sinh mh li^ ? x. ^ 


Bo- = -f 


B,. = _ B, . 

3i 2P li 
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4i 2U sinh mh li 


3p 

{ (cosh nili-l) + h } 


(8 .13) 

(6 .14) 


where 


1 

. X 2 p 

i = 1,2, x^ = X, X2 = z, II = (-“y 4S) J m j , a 


(~f 


( b .15) 


(iii) Generalized Reynolds Biuation 

Integrating the equation of continuity 

(8 . 16 ) 

across the film, we get the generalized Reynolds equation as 


3u 3v _3w 
3 X 3 y 3 z 


3p 


3p 


3x ^ 3z { P^^(ll,il,h, 4^^) 


“ ir' ^ ^ a^2^11,il,h,4^)} + V + 7 


(8 .17) 


where 7 is the relative normal velocity of the surfaces and 7^ is 
the velocity of the fluid towards the porous region at the cartilage 


surface and 
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i = 1,2 
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^ ^ /'ET G \ h^ , , 2, 11 S, h^ NH 

ffl = + £ h - coth ^ 

(8 .20) 

^21 - = . t 
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(8 .21) 

^31 "" 11 
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How the prohlem remains to find the Telocity of the 
fluid at the porous surface towards the porous region* For this 
purpose, let us recall that the hyaluronic acid i s a long chain 
polymer compound. The molecular weight of hyaluronic acid molecule 
is 5 ** 10^ (aijproximately) and molecular length of the order 
5 X lO"^ - 10"^ cm, Bowson [115 ] • These molecules normallj*- 
can not pass throu^ the porous matrix which is the basis of 
boosted lubrication. Thus it can be assumed that the flow 
through the cartilage matrix is Newtonian, and hence the velocity 
of the lubricant throu^ the porous matrix shall be governed 
by the usual Darcy’s law [ 135 1 • 


Thus [ 135 3 » 
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Y 
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6 

a I 

M 3 y 'y=h 


(b ,23) 


Prom the requirement of continuity we have for the porous matrix 

^ o 


V.q = - 


p = o 


(8 .24) 


so that, since 



V 


2 

P = 


o 


(8.25) 


Ihe problem now ireduces to get the solution of the Reynolds 

equation (s .1?) for the pressure in the oil film, simultaneously 

with that of Iiaplace eqn, ( 8 . 25 ), for the porous matrix, with a 
3P 

common -r- at the boundary, 
cty 

J^ssuming the thickness of porous matrix to be very small 
as coE^iared to the radius of the sphere R) , the average 

pressure at any radial section of the porous natrix can be taken 
equal to the pressure in the lubricating film. It is further 
assumed that is linear across the matiix and is zero at the 

outer surface of the porous bearing shell. 


This assumption leads to 
2 

— ^ = constant = K^(say) 

3y 


(8.26) 


H^ce from Laplace equation 
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Integrating (8*26) j 'we get 


op 

37 


= K^y + ^2 


(8 . 28 ) 


where Kg constant of integration. 

0 

With the boundary condition ^ = o at y = h + we get 


3P 

9y 


= K, (y-h~H^) 


( 3 . 29 ) 


Hence 


OP 


9 7ly=h 




Substituting the value of K^ from eq.n. ( 8 . 27 ), we get 


OP 
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, 0 o^- 

= H^( -f + — 

y=h Ox Oz 


(8 .50) 


(8 . 51 ) 


Thus 


4> H 2 ^2 

P ^ Ox^ Oz 


Substituting the value of in e<in. (8.1?) and simplifying, we 
get 


( 8 . 52 ) 
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Equation (8.55) is the genoralized Rejmolds e<lua,tion governing 
the pressure distribution in the cavity of the joint. 


8.4 NOBITiiL MO^riQg OF THE Em;IAi: JOUTO 

Attention is norf restricted to the stady of norisal 
motion of the joints, i’herofore, the effects arising out of 
sliding are neglected, i.e. == = o, i = 1,2, Thus eq.n. 

(8 . 55 ) takes the form 




[ { + 
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(8.54) 


The film thickness is 


h = c (l-e cos 6 ) 


and the velocity of squeezing is 
dh dE 

Y = ~ = -o ■rr— cos 6 
dt dt 


Using the spherical polar co-ordinates for the system under 
considomtion, i.e, using 

X = R sin 6 cos C 


y = R cos C 
z = R sin 9 sin C 
eqn. (8.54) assumes the form 


(8.55) 


(s .36) 


( 8 .37) 

(8 .38) 

(8.39) 
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11 

de 


dt 


cos 0 


( 8 .40) 


It is reasonable to assnmo = o (i.o. the pressure is consta-nt 

d Z 

along Z direction). V/ith this assunption, eq.n, (8 » 40 ) reduces 
to 

— [ sin e{F,, (N,£,h,<)J + — 


t,2 . . 90 

R sin 0 


' 11 ^ ’ * *^1 
d£ 

= - c rT“ oos 6 
dt 


(8A1) 


To solve the Reynolds eqn. (8»4l) the pressure boundary 
conditicffi must be known. Since the pressure reaches a maximum 
at 6 = o, therefore the first condition becomes 


= 


d0 


o at 6 


(8 . 42 ) 


However the other condition is debatable as the exact 
extension of lower cavity of bone is undetermined. If one asaimes 

the lower cavity to termimte at some points = + S ^ then the second 
condition would become 

p s= o at 0 — 4- 6 (s # 45 ) 

But this again involves an unknown parameter 6 , However it is 
reasonable to assume that ® ~ (approximately) . This assun^ttion 
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may be justified because the pressure reaches a peak at 0 = o and 
falls very sharply on both sides. Hence the influcaace of the 
boundary on pressure distribution is likely to be negligible 
after some distance from 0= o. 

Thus we have the second boundari,^ condition as 

p = oat 6 = +-2 (8.44) 

Integrating eq.n, (8.41) tasing the boundary conditions 
(8.42) and (8 44) the equations governing the pressure gradient 
and pressure distribution! are, respectively 

. -|li E^c U sine X 

, < f,;,) 

-- hh + ^ 31 ’ * ^31 ‘hk ( i- - 
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The exjsression for load capacity is 
2 

W = 2'n' R j p sin 6 cos 0 de 
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(8 .47) 


or 
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Ron-dim ensionalizing, rising 
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where and 4 are characterization of slip and porosity effects, 
p s 

respectively. Thus the non-dimensional load capacity is 
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where, 
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The time of approach, t, from an initial film thickness to a 
final film thickness, in terms of corresponding eccentricity 
ratios, is given hy 

t- rlli. 
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The non-dimensional time of approach is 
2c^t 1 
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8#5 BIBULTS nSCTJSSIOII 

(i) Djiaens ionless Parameters 

Apart from the usual dimensionless numbers encountered in 

the classical theory four new parameters d > d » ^ and L, are of 

s p 

importance in the study of the present joint lubrication problem, 

fhe slip parauBter d the permeability d have since long 

s p 

been identified and discussed. It would suffice here to maition 

that a larger value of d 02: d signifies a greater slip or a 

s p 

greater permeability and a value of dg ~ ° '!*p " ° signifies 

no slip or a non-porous body. It may be noted that for a non-porous 
body the slip is necessarily zero, however, a zero slip does not 
imply a non-porous case, 

The parameters W and L have already been discussed in 
Chapter III, In the present context, it suffices to say that, 1 
and L represent concentration and length parameter of the hyaluronic 
acid molecules, respectively, in synovial fluid. Higher is the 
value of N hi^er is the concentration of the molecules, and lower 
is the value of L higher is the chain length of molecules, Thus 
the joint lubrication characteristics would be prominent when N 
is higher and L is lower, 

(ii) Diseased Joints 

Ihe properties of the synovial fluid undergo a marked variation 
with age, specially in case of a diseased or damaged joint. In a 
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rheumatoid arthritic patient, the s;;uiovial fluid xna;.’- even Icose 
its non-Newtonian behaviour, Bloch and Dintenfass [ 134 ] • This is 
characterized by the cases of II = o or L -+ o= , These cases actually 
correspond to the classical theory approach . Burch et . al [135 j 
indicated that the concentration of hyaluronic acid molecule in 
synovial fluid decreases for a diseased joint. For old joints, the 
cartilage becomes softer, less resilient and may even vfear out in 
certain cases. These two cases can be diaracterized by a decreasing 
value of N and d respectively. Furthermore, it is seen that for a 
rheumatoid joint hyaluronic acid molecule has a lower molecular 
weight » A lower molecular wei^t corresponds to smaller length 
of the molecule. This, in this study, is depicted by a larger value 
of L. 

Thus it is evident that, the present stut^ is capable, 

throu^ the parameters 4 > 4 » ^ and N , to analyse the 

s p 

behaviour of a diseased joint also, 

(iii) Joint Characteristics 

The non-dimensiaxal load capacity t is given by equation (8.50) • 
Even a casual glance, tells that the slip parameter 4^ is multiplied 
by the factor (l-H^), in the numerator as well as the denominator. 

The same is true for the expression of the dimensionless time of 
approach!, eqn. (8.55). This is an important result. This indicates 
that as N increases the effect of the slip will decrease, because 
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2 

*5*2 (l-l^ ) will deci'oase, Incraaso in IT, is a conseiuence of the 
Boosted lubrication. So as ths concentration of th: hyaluronic 
acid molecules increases, the effective viscosity bocories more and 
the effect of slip on the load and response time decreases. 

Curves for dimensionless load capacity T? vs. the various 
parameters involved are shown in Figures 8 .3 to 8 .9 . Figures (s .5) 
to (8.5 ) s33DW that the load decreases as the value of slip parameter 
increases. The qualitative trend is similar to that obtained by 
the classical theory approach. Of importance is the fact that load 
decreases as (a) the porosity increases j Fig, 8,3, (h) the concentration 
decreases ; Fig, 8.4 and (c) the chain length of hyaluronic acid 
molecule decreases | Fig. 8*5* These results establish theoretically 
that a diseased or a damaged joint will carry a lesser load. 


Graphs for the load vs . (}> ; Figs . 8 .6 - 8 .7, indicate 

]P 

a similar trend that the load decreases as 6 increases for a 

P 

fixed value of other parameters. Fig, 8,7 5 which is a plot of 
W vs, for different values of H, indicates that if 6 
increases and H decreases, simultaneously, the decrease in load 
will be veiy rapid. 


Fig, 8,8 is a graph of ^ vs, H for the porous (indicated 
by continuous lines) and the non-porous (dotted lines) cases for 
various values of L, The curve for L-*- <» gives the minimum load 
in both cases. This case corresponds either to the classical 
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theo3?y approach or the reverse case of diseased ioint ’."here the 
synovial fluid might have lost its non-Hewtonian characteristics. 
Obviously this curve is straight line, since it no longer is a 
function of II, Also, since the porosity of the cartilage has been 
established beyond doubt, the case of non-porous cuirves, in spite 
of their higher load capacity, are not of much significance. These 
curves are however given here only to illustrate the significant 
q.uantitatiTe difference in the load capacity. 

Dimensionless load W vs, L(Pig. 8*9) shows a very slow 
decrease in the load for low values of N, Indicating that for 
low concentration of the hyaluronic acid molecules in the synovial 
fluid the load is very low, even when the chain length of the 
molecule may be high. But as the concentration further increases, 
the capacity becomes higher for long chain length (low value of l) 
of the molecule. However, as the chain length decreases, the load 
capacity decreases rapidly, irrespective of the value of the 
concentration. This happais in the dic?eased rheumatoid joint. 

The curves for the film thickness vs , the time of approach 
are shown in Pigs, 8, 10-8. 13 • The study of these curves is very 
vital for the joint lubricatim problan because this would indicate 
the time it would take the two cartilage surfaces to coae into 
intinate contact, which actually will happen if the film thickness 
become zero, i*e* all the fluid is squeezed out from between the 


two surfaces 
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It is seen that the tine of approach is sufficiently large 

for low values of the slip parameter ; Pig, 8.10 , hut an increase 

in the slip decreases the time of approach ; thus preserving the 

lualitativB trend of the classical theory approach, A similar 

result is visible from Pig. 8,11 ; that I decreases as the porosity 

increases. The time of approach for the slip case (continuous lines) 

and the no-slip case (dotted lines) is almost identical. The 

most important aspect of this figure is the fact that as the value 

of increases a finite time of approach is obtained. This becomes 

sufficiently small for high valie of ({j , which indicates that all 

P 

of the synovial fluid is either squeezed out or passes into the 
pores, causing an intimate contact of the cartilage surfaces in a 
very short period of time. This may be a cause of consistent severe 
pain in the diseased rheumatoid joints. 

Almost similar results are obvious from Pigs . 8 .12 - 8 ,15 » 
wherein the results for the non-porous case are indicated by dotted 
lines and for that of the porous case by continuous lines. The time 
of approach is seen to decrease as either the chain length of the 
hyaluronic acid molecule decreases or as the concoatration decreases 
(i,e. as L increases and N decreases) . Finite value of T are once 
again indicated for high value of L (low ciain length) and low value 
of IT (low concentration) , 
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time of aoproach for different values of coupling number 


207 


GELETER IX 


SLIIIftJiy jJ-TD C0ITCLUSI02IS 


The classical lubrication theory is haseQ. upon a lundametital 
assumption that the fluid film behaves in a lTe’.vtonian '«vay, i»c. 
on the classical continuum theory, Hoiirever, the atoaistic models 
of the material, especially of fluids that possess coarse structure 
and fibres (such as colloidal fluids, liquid crystals and fluid 
containing additives) have shown that the mass density can fluctuate 
violently with the size of the volume element, when it is below 
a certain critical limit but still in the macroscopic range. Thua 
the continuum assumption for mass density is no longer applicable. 

Another important class of problems where the classical 
theory is likely to give erroneous results is when the characteristic 
length scale is comparable to the average grain or molecule size contained 
in the medium. In such situations the molecular or the granular 
constituents of the medium are excited individually and the 
intrinsic motion of the material constituents must be taken into 
account, Typical of this class are the problem of flow through 
narrow passages, such as those found in human arteries and in 
lubrication theory. 

The classical lubrication theory also contains no mechanism 
to explain the presence of the residual film so ccmmonly observed 
by various workers [ 22,28,29,51,56 ] , Thus there is a strong 



208 


motivation for extending the range of e-pplicahility of continunm 
mechanica, ya—f-icularjj in the s'h.idy of fluid fio^r of the tj-pe laaationed 
abovejusing tne field theory approach to describe the nacrcocopio 
manifestation of nicroscopic events, i,e, nicroactions and deform aticns. 

Kicroscopic effects generated by the nicroaotions of particles 
in suspension of viscous fluids, or bjr the uxolecules in the proximity 
of a solid surface, drastically change the character of flew between 
narrow passages. The origin of the rheological abnormalities is then 
to be investigated on the basis of the theory of micro-mechanical 
motions . 

One way of accounting for these micro-mechanical motions is 
through the so called theory of fluid miorocontinua, advanced by 
Eringen [ 40} » which the continuous media are now regarded as 

sets of structured particles, which not only contain mass and 
velocity but also a substructure, i.e. each material volume element 
contains microvolume elements which can translate, rota+e and 
defoirm independently of the motion of the macrovolume. Thus a 
mechanism is provided in this theory to treat materials which are 
capable of supportir^ local stress moments and bo^ monents ani in 
addition are influ aaced by the mioroelanent spin inertia. 

The micro fluid theory is however too complicated and the underlying 
theoretical problem is not easily amenable to the solution of non-trivial 
problems in this field. This led Eringen to postulate a subclass of 
these microfluids, called micropolar fluids [ 44 ] » which would still 
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exhi'bit effects arising from particle micronotion. For this class 
of fluids, the defornBtion of the fluid microelement is ignored, 
nevertheless, microrotation effects are still present and surface 
and body couples are permitted. 

Physically the theory of micropolar fluids may serve as a 
satisfactory model for description of the flow "behaviour of polymeric 
fluids and fluid suspensions, or wherever the microrotation my 
play an important role. One is then naturally inclined to apply 
this concept of the microcontinuum to'^^i© analysis of a wide va-riety 
of fluid flow pro"blems wherever a further refined model of fluid 
"behaviour is sou^t . 

Motivated "by these considerations tte lubrication theory for 
micropolar fluids has been developed by many authors in the last decade, 
"Various traditional lubrication problems have been reported. The 
results obtained are encouraging and illustrate the capabilities of 
this theory, in explaining the phenomenon of the eiiliancement of 
effective viscosity in thin films. Even after a decade, this study 
seems to be in a preliminary state. Scores of lubrication problems 
have yet to be solved from this view point • For example, little work 
has appeared on dynamically loaded bearings and porous bearings from 
the microocntinuum view point, Iforeover, the problems of finite 
dimensional bearings and bearings with rough surfaces from this view- 
point are yet to be analyzed. 
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Htus, ir. this thesis, these luhrication probleins are studied froc 
the microccatin.u'u^i point of vie'.v, in an atter.pt to erplain theoretically 
the rheolcEpical alnorralities so frecpuently observed in fluids v/ith 
additiveo, especially v/hen confined to narrow passages. 

The thesis begins with an introductory type Chapter, eaphasizing 
the existence of the residual filn in the proximity of a solid surface 
and laying the foundation for the application of the theory of micropolar 
fluids to certain hitherto untouched lubrication problems. Hien in 
Chapter II, a generalized Reynolds equation applicable to finite lubrication 
problems, is derived, using the micropolar fluid theory. The equations 
obtained in Chapter II are applied in Chapter III, to a two-dimensional 
problem of squeeze film of a ball in a spherical seat and to souk three 
dimensional non-cyclic squeeze films, assuming the characteristic 
coefficients to be constant, in an attempt to study the effects of rigid 
particle additives for the tnree dimensional micropolarity model. 

Increase in effective viscosity due to the micropolarlty is established 
theoretically. It is also shown that the theoretical effects of the 
micro polarity on three dimensional lubrication are identical to the 
two dimensional problems, at least qualitatively. 

A very important result of this study is the dependence of 

bearing characteristics on a non-dimensional length parameter 

0 *V X /2 

l(= j , y/here c is the bearing clearance and £ = ( , a fluid 

property having the dimensions of length) . This parameter links the 
fluid with the bearing geomd:ry and can be thou^t of as a number which 
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is inversely proportional oithor to the length of the nolecnle of 3, 
non-lTevrconian luhricant or to the inoleoular sire of ohe addiilve. 

It is foau-l thafc snaller the L(i^. larger the , or snaller the 
clearance) more prononnced are the micropolar effects. 

“his dependence of hearir^ characteristics, especially effective 
viscosity, on the geometry of the system, may also give a possible 
explanation for discrepancies in actual viscosity measurements made 
using viscometers of different designs or sizes. 

It should also he noted that most of the earlier workers have 
2 

used the relation o 1. II <1. However, in Chapter III it is shown that 
o 5. results of those workers are grossly exag’gerated 
due to non-admissihle values of N. 

The next problem (Chapter IV), deals with the dynamically loaded 
short journal hearings in which’ both wedge and sq.ueeze films are 
active. The Reynolds equation for the general case of dynamically 
loaded infinitely short hearing is derived, where the lubricant 
is assumed to he micropolar* Detailed consideration is given to 
the dynamic hdiaviour of squeeze film in a short journal hearing 
under a sinusoidal load with no journal rotation. Various hearing 
characteristics are obtained assuming a full film to exist, Hie 
micropolarity of the fluid results in more i^sistance to journal 
motion, thereby allovFing for smaller eccentricities for a constant 
load. The overall oonclusion of this study is an increase in effective 
viscosity due to the micropolarity of the lubricant , 



2.12 


The next protlem (Chaptei’ 7), is tliat of the roller hearings 
ill comhined rolling, sliding and normal motion analysed u'lder 
cavitation houndar;;.' conditions. Various hearing characteristics 
are obtained and represented grrphically . The overall conclusion 
of the study, in confomity y,'itla the erperimentc.! results, is an 
increase in effective viscosity, resulting' in an increase in load 
carrying capacity and a decrease in coefficient of friction, when 
analysed from microcontinuum view point. 

It has been experimentally observed [ 27 ^ that there is a 
sudden decrease in the coefficient of friction on addition of 
long-chain polar compounds . A graph of the values of coefficient 
of friction plotted against molecular weight (see ref. 156 » P*550) 

shows the greater the number of atoass and loiter the molecule the 
lower the coefficient of friction will be. Holding h constant 

h 

(whei^e is the minimum film thickness) , and recalling that L = 

(L is defined in this way in Chapter V), L for the curves in Figs ,5. if 
to is seen therefore to depend only on the fluid, that is on 

the i material characteristic length. Considering the length of 
molecules of the boundaiy lubricant to be proportional to the material 
characteristic length, it is observed that the portion of the curve 
(n = 0,7) in Figs, 5,^ to 5,8, to the right of the minimum satisfies 
the result that an increase in the length of molecule (coCTes ponding '' 
to a decrease in L) causes a decrease in coefficient of friction and 


V ice-versa 
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In ■':he next pro’olen (Chapter '"l) j micropolar fluid tlieorj' is 
ap’-lied to study the effects of surface roughness. Generalised 
forns of Le^molds eiuation are derived for three cases of roughness, 
iie-nely, transverse rcaghness, longitudinal roughness and isotroric 
roughjaess, using, t'ne stooiiastic approach. For mtheniatical sinrlicity 
only the stationery part of the hearing' is assuaed to he rotigh. This 
theory is suhseq,uent 2y applied to the prohlem of infinitely long 
journal hearing in order to study the interaction of micro polarity 
with surface asperities, using the half Sommerfeld boundary conditions. 
Various hearing characteristics are obtained nunerically and results 
are compared \'dLth smooth hearing using nominal film thiclmess concept. 
The results obtained are in qualitative agreement with ilewtonian 
results for rough surfaces. However, it is shown that the effects 
of surface roughness are more pronounced, when analyzed from 
microcontinuum view point • 

A number of experiments in which the evidence of rheological 
abnormalities, liJce enha-ncement of viscosity, has been rather 
extensive and convincing, are of squeeze film type [ 22,28, 29f31j 563 • 

It has been suggested that the asperities on the solid 
surface might account for the rheological abnormalities and iresidual 
film [ 50 ] . 

Motivated by these considerations, the next problem (Chapter VIl) 
deals with the squeezing between rough parallel plates. Three standard 
geometries, namely, infinitely long parallel plates, circular plates 
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and finite rlateS} are considered. Varici^s r-neene xilr. 

characteristics arc- obtained for the t'nree esses *f r-'njrh.necs . The 
results are corrared rith the correspendirr* smooth plates. The 
rather important peonetrj, i.e. the finite rectan2"lar vla'^es saueeze 
film is analysed in details. Since this geometry can distinguish 
betv.'een the one-dimensional rougl'niess and isotropic rougliness, a 
comparision is made between these tn’o cases. Further, for a constant 
area bearing, it is observed that the aspect ratio (width/length) 
plays a vital role for one-dimensional roughness, i.e. the bearing 
characteristics, for a given roughness, may increase or decrease 
depending upon this parameter. 

The last problem (Chapter Vlll) is of synovial joint 
lubrication studied from microcontinuum view point. The syaovial 
joints provided by the nature in the human body to carry out the 
trouble free motion of one bone past another, have long been 
identified as bearing systems. The present situation regarding 
joint lubrication is, however, not very clear. Earlier werk, 
based upon the class ioal continuum theory, could not give a 
satisfactory explanation for the increased effective viscosity 
in synovial joints. The miorocontinuum approach represents a more 
realistic model for explaining the conplex mechanism occuring in 
human joints. 

The bdmviour of the synovial fluid, which has long chain 
hyaluronic acid molecules, is oensidered to be governed by the 
micro polar fluid theory. The porosity of the cartilage surface 
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and the slip at the porous boundary are taben into account and a 
modified form of Heynolds equation is derived. i?he joints are 
approximated by a spherical bearing and various bearing diaracteristics 
are obtained in d-.-tails . The salient feature of Lhis study, is an 
increase in effective viscosity. This theory also establislies 
theoretically, that for a diseased or a damaged joint the load 
capacity ana the response time decrease • This decrease may result 
in pain in ca-se of diseased rheumatoid joint , 

The correlation of the results obtained in this thesis 
with various experimental studies may perhaps be believed coincidental 
and is definitely open to criticism since, in the light of these studies 
no experimeatal verification of the micropolar theory has been made. 
Nevertheless, at this early stage such criticisms are natural and in 
fact healtliy for a crucial scrutiny of the theory and its ultimate 
acceptance. The results obtained are, however, interesting, and 
it is our strong belief that the micropolar fluid theory has the 
potentialities to explain the rheological anomalies in fluids 
containing additives through the parameters N and L when they flow 
through narrow passages, 

Although the various results based on this theory are in 
perfect qualitative agreement with tlie experimentally observed results, 
an attsnpt for a quantitative fit would not be of any significance 
(since two undetermined paramcfeers exist for curve fitting) unless some 
numerical values could be ascribed to the coefficient X and Y , 
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A sigaifioant step in this direction is hy Arirr-n et.al [ 137 j 
who estahlished a relationship bet. Teen the nlc_cpoiar oocificients 
and henatocrit concentration in blood sanples. Further, Inpliarello 
and Sevilla [ 138 j , GoldsEith and lH&son [ 139 ] » Sak-haveva 
et. al. [ 140] and ?uks [ 28,29 ] estatlisrai a reLationshii bot'Teen 
the effective viscosity and the molecular 'weight (particle size) • 

It is hoped that these studies would ra-ise the curiosity 
of serious research workers to this rich field. The most useful 
contribution in this direction would be a repetition of the past 
experimental st'udies, in the light of these microcontinuum effects 
and thereby providing a much needed correlation with the theory. 
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